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The non-Abelian topological order for superconductors is characterized by the existence of zero- 
energy Majorana fermions in edges of systems and in a vortex of a macroscopic condensate, which 
' obey the non-Abelian statistics. This paper is devoted to an extensive study on the non-Abelian 

, topological phase of spin-singlet superconductors with the Rashba spin-orbit interaction proposed 

in our previous letter [M. Sato, Y. Takahashi, and S. Fujimoto, Phys. Rev. Lett. 103, 020401 
(jj ' (2009)]. We mainly consider the s-wave pairing state and the d -\- id pairing state. In the case of 

\ d -\- id-wave pairing, Majorana fermions appear in almost all parameter regions of the mixed state 

■ under an applied magnetic field, provided that the Fermi level crosses fc-points in the vicinity of the 

r point or the M point in the Brillouin zone, while in the case of s-wave pairing, a strong magnetic 
field, the Zeeman energy of which is larger than the superconducting gap is required to realize the 
(-H ' , topological phase. We clarify that Majorana fermions in Rashba spin-singlet superconductors are 

Q ' much more stable than those realized in spin-triplet p + ip superconductors in certain parameter 

^ ' regions. We also investigate the topological number which ensures the topological stability of the 

phase in detail. Furthermore, as a byproduct, we found that topological order is also realized in 
conventional spin (or charge) density wave states with the Rashba spin-orbit interaction, for which 
' massless Dirac fermions appear in the edge of the systems and charge fractionalization occurs. 
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I. INTRODUCTION 

Q ■ Topological states of condensed matter systems are characterized by a bulk topological number such as the Chern 
Q ' number (or the TKNN number) which represents a topologically non-trivial structure of the many-body Hilbert space. 
In such phases, topologically protected surface states and fractionalized quasiparticles, e.g. anyons, appear.l"''^l jj^ 
particular, when topological order is realized in a certain class of superconductors, this topological phase supports 
^ : the existence of chiral Majorana edge modes and a Majorana fermion in a vortex core.—"— Vortices with Majorana 
' fermion modes are neither fermions nor bosons, but non-Abelian anyons, obeying the non-Abelian statistics for which 
00 . the exchange operations of particles are not commutative;3.".2. Because of this remarkable feature, a vortex with a 
' Majorana fermion may be utilized as a decoherence-free qubit, and plays an important role for the realization of fault- 
' tolerant topological quantum computation.—"— The state with non-Abelian anyons, which is called the non-Abelian 
topological phase, has been discussed to be realized in the fractional quantum Hall effect state with v = 5/2 and 
[ 12/5<^"— It has been also known that spin-triplet superconductors such as chiral p + ip superconductors ^i^^"— 
. and noncentrosymmetric (NCS) p-wave superconductors with broken time-reversal symmetry^^i possess a zero energy 

■ Majorana mode, and realizes a non-Abelian topological phase. In general, fully-gapped spin-triplet superconductors 
^ , support non-Abelian anyons if the number of the connected Fermi surfaces are odd, in the case without time-reversal 

symmetry.^^ For the spin-singlet s-wave superconducting state, it was pointed out by Fu and Kane that non-Abelian 
^ . , anyons are realized in the proximity with a topological insulator Also, the non-Abelian anyons in the s-wave pairing 

■ state was discussed before in the context of Axion strings in cosmological systems. ^'''^^ 

Recently, the present authors proposed another scenario of a non-Abelian topological phase in NCS s-wave super- 
fluids or superconductors. We pointed out that in the presence of the Rashba spin-orbit (SO) interaction, s-wave 
superconducting states show a transition to the non-Abelian topological phase with non-zero Chern number, under 
an applied strong Zeeman magnetic field^SS,. Also, independently, it was proposed by Sau et al. that such systems 
can be realized in heterostructure semiconductor devicesi^ The idea was subsequently generalized by Aliceai^ 

In this paper, we explore extensively properties of the non-Abelian topological phase realized in NCS spin-singlet 
superconductors with the Rashba SO interaction, which was considered in our previous letter There are two main 
purposes. The first one is to present the detail analysis of chiral Majorana edge states and a Majorana fermion 
mode in a vortex core in the case of the NCS s-wave superconductor, and the calculation of the topological number, 
which are omitted in ref.^®. The topological order for time-reversal symmetry broken systems in two dimensions is 
characterized by the first Chern number. We present a formulation for the calculation of the Chern number with the 
use of a winding number, which makes the estimation of the topological number easier. Using this formulation, we 
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explore the non-Abelian topological order realized in NCS spin-singlet superconductors. We, furthermore, analyze 
the vortex core state by solving the Bogoliubov-de-Gennes (BdG) equation. We obtain the zero energy Majorana 
fermion solution when the Zeeman energy due to an applied magnetic field is larger than the superconducting gap 
A. We also discuss that, in some parameter regions, the Majorana fermion in NCS spin-singlet superconductors is 
remarkably stable compared to that in chiral p + ip superconductors. This feature is crucially important for the 
application to the topological quantum computation. The superior stability of Majorana mode in a vortex core of 
NCS spin-singlet superconductors stems from the fact that when the Zeeman energy p-bHz satisfies the condition 
A > /iB^z — A > 0, the Majorana fermion is mainly formed by the superposition of quasiparticles in the vicinity of 
the r point (or the M point) with the Fermi momentum kp ^ (or (tt, tt)), in the long-distance asymptotic regime 
sufficiently far from the center of the vortex core. Because of this property, the excitation energy Eq in a vortex 
core which separates the zero energy Majorana mode and the first excited state is much larger than a typical energy 
scale of the Andreev bound state of vortex cores ^ A'^/Ep. Furthermore, the vanishing Fermi momentum for a 
Majorana fermion implies that decoherence due to quantum oscillations of quasiparticle energy with a period ~ 1/kp 
raised by inter-vortex tunneling,— which may be an obstruction for the implementation of the topological quantum 
computation, is substantially suppressed. 

The second purpose of this paper is to extend the scenario of the non-Abelian topological order for the case of 
s-wave pairing state to other spin-singlet pairing states. In particular, we consider the cases of d + id-wave pairing, 
for which there is a full-gap in the energy spectrum, ensuring the nonzero Chern number. It is demonstrated that in 
the d + id-wave pairing state with the Rashba SO interaction, when the Fermi level crosses /c-points in the vicinity 
of the F point or the M point in the Brillouin zone, the non-Abelian topological order, which supports the existence 
of chiral Majorana edge states and a Majorana fermion mode in a vortex core, appears under an applied magnetic 
field. In contrast to the case of s-wave pairing considered in refi^^— , for which Majorana fermions appear only when 
there is Zeeman splitting larger than the superconducting gap, a small magnetic field larger than the lower critical 
field suffices for the realization of the non-Abelian topological order in the d + id-wave pairing state. Thus, it may be 
easier to realize Majorana fermions in the NCS d -|- id-wave superconductor than in the NCS s-wave superconductor. 

Furthermore, we consider another direction of the extension of the scenario for the non-Abelian topological order. 
Our results for Rashba s-wave superconductors imply that the topological order is also realizable in the conventional 
spin density wave (SDW) state or the charge density wave (CDW) state with the Rashba SO interaction. We 
demonstrate that in these density wave states, the Abelian topological order appear under applied magnetic fields, 
leading to the existence of gapless edge states described by the Dirac fermion, which is analogous to the surface states 
of the topological insulator. We also discuss the scenario of charge fractionalization in the topological density wave 
states. 

The organization of this paper is as follows. From Sec|n]to Sec lIIIi we introduce the model for superconductors 
with the Rashba SO interaction in two dimensions, upon which our analysis is focused, and, as a first step of our 
analysis, classify the parameter regions of the model, in which different topological phases may be realized. In Sec lIVi 
we explain the duality relation which holds for our model Hamiltonian. This duality relation was utilized for the 
analysis of topological properties in ref.-^^. In the most part of this paper, we do not use the duality relation, but 
instead, confirm the argument based on it developed in ref.^^ by adopting a more direct approach to this issue. In 
Sec|Vl we analyze and discuss the topological number characterizing the non-Abelian topological phases realized in 
NCS spin-singlet superconductors. In particular, we prove the relation between the Chern number and the winding 
number, which is useful for the investigation of the topological order. Using this relation, we elucidate the general 
condition for the realization of the non-Abelian topological order. On the basis of the analysis of the topological 
number, we obtain the phase diagram of the topological order for spin-singlet NCS superconductors. We also give 
some physical arguments on the origin of the topological order in spin-singlet NCS superconductors. In Sec lVIl the 
numerical results for chiral Majorana edge modes are presented for the cases of s-wave pairing and d -I- id- wave pairing. 
In Sec lVIIi we consider an approximated but analytical solution of the BdG equation for Majorana zero energy modes 
in vortex cores. We, also, discuss the superior stability of the Majorana mode in NCS spin-singlet superconductors 
compared to that in chiral p + ip superconductors. In Sec lVIlH topological density wave states in which gapless 
Dirac fermions on the edge of systems appear and charge fractionalization occurs are considered. In Sec lIXl we give 
a summary of our results, and also discuss possible realization of the NCS spin-singlet superconductors with the 
non-Abelian topological order in real systems. 

Some technical details are presented in Appendices. In Appendices [X] and |B1 we derive useful formulas for the 
Chern numbers and the winding numbers, which were used for the discussion on the topological number in SecfVl 
Supplementary discussions related to the topological argument given in Sec|V] are presented in Appendix [C] The 
details of the derivation of the BdG equation for a singlet vortex are given in Appendix |DJ In Appendix |E1 we discuss 
another mechanism of non-Abelian anyons in spin-singlet superconductors, which was first discussed in ref.^"*, where 
non-Abelian anyons are realized in time-reversal invariant s-wave superconducting state without a Zeeman magnetic 
field. This discussion is relevant to the non-Abelian topological order realized in an interface between an s-wave 
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superconductor and a time-reversal invariant topological insulator proposed by Fu and Kanei^ 

II. MODEL 

In this paper, we consider spin-singlet superconductors with the Rashba SO interaction'^*' in two dimensions. For 
concreteness, we define our model in the square lattice, while the following argument does not rely on the particular 
choice of the crystal structure. It is also noted that our analysis and results are also generalized straightforwardly to 
other type of anti-symmetric SO interactions raised by the lack of inversion center of systems. The Hamiltonian is 
given by 

,0" k,a.(7' k,(T,(T' 

+ ^ E ^--'WcLcLfc,, +i E K'Ak)c^kaCka', (1) 
k.{T,a' k,a.{T' 

where {cka) is a creation (an annihilation) operator for an electron with momentum k — (kx,ky), spin a. The 
energy band dispersion is e(fc) — —2t{coskx + cosfcy) — /i with the hopping parameter t and the chemical potential 
/i, and the Rashba SO coupling is aCo{k) — a{sin ky, — sin k^) {a > 0). We also introduce the Zeeman coupling 

-fJ'BHzJ2k,aa'('^^)'^'^'4:a^ka' in the Hamiltoniau. 

For the spin-singlet superconductors, the gap function A^a' (k) is written as 

A,,,{k) ^iA{k){ay)aa, (2) 

with the y-component of the Pauli matrices ctj {i = x, y, z). In the following, we assume two different full-gapped spin- 
singlet superconductors: the first one is the s-wave pairing, A{k) = Ag, and the other is the d + id-wave pairing. For 
the zd-pairing, we consider two possible realization on the lattice, A(fc) = A^^^ (cos ky — cos k^) + iA''^^ sin k^ sin ky, 

or A(fc) = A^^"* (sin^ fcj: — sin^ ky) + zA^^-* sinfca; sinfcj,. The amplitudes A^ and A^j*"* (i = 1,2) are chosen as real 
and positive. The second type of the d + id-wave pairing includes higher harmonic contributions, which may arise 
depending on detailed structures of electronic bands and the pairing interactions. We use these two types of the d+id- 
wave gap to clarify that, although both of them support the non-Abelian topological order, the precise condition for 
the non-Abelian phase slightly depends on the detail of the gap structure. 

For a noncentrosymmetric superconductor, the parity mixing of the gap function generally occurs^ir— Therefore, 
in addition to the spin-singlet component of the gap function, the spin-triplet one is induced generally. However, if 
the spin-singlet amplitude dominates the gap function, the topological nature is not affected by the spin-triplet one. 
We neglect the spin-triplet component in the following. 

In the following, we mainly consider the case that SO interaction is much larger than the Zeeman energy; i.e. 
a|£o(fc)| ^ (J-bHz, which is an important condition for the stability of the superconducting state against the Pauli 
depairing effect due to the magnetic fields. 

III. GAP CLOSING CONDITION 

In general, continuous topological phase transitions between topologically distinct phases occur only when the 
energy gap of the bulk spectrum closes. Thus, to identify parameter regions for which different topological phases are 
realized, we first examine the bulk spectrum of the system. To obtain the bulk spectrum, we rewrite the Hamiltonian 
as 



(3) 



where the BdG Hamiltonian 'H(fc) is given by 



%l(h\-\ s{k) ~ i^BHzCFz + aCo{k) ■ a iA{k)ay > , , 

"■^''>-\ -iA{k)*ay -e{k) + fisHziz + a£o{k) ■ cr* I ■ 



Diagonalizing the BdG Hamiltonian, we find 



E{k) = ^e{ky + a^Coiky + n^H^ + |A(fc)|2 ± 2^eikya^Co{k)^ + {e{kY + | A(fc)|2)^B-ff|. (5) 
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In our model, the gap of the system closes only when the following condition is satisfied, 

e(fc)2 + a'Coikf + ^1-^' + |A(fc)|2 = 2^ eikya^Co{ky + (e(fc)2 + |A(fc)|2)^2 (g) 

From a straightforward calculatioit^, it is found that this condition is equivalent to 

e{kf + \A{k)\^=filH^^+a^£„{kf, \A{k)\a£„{k) = 0. (7) 

We examine the gap closing condition using ([7]) for the s-wave pairing state and the d + id-wave pairing state in the 
following. 

A. s-wave Rashba superconductor 

For the s-wave pairing, the second equation in (O is met only when Co{k) = 0. Therefore, the gap closes at 
k = (0, 0), (0, tt), (tt, 0), (tt, tt). Substituting those values into the first equation in we have three different gap 
closing conditions. ( In the square lattice, the condition at fc = (tt, 0) and that at fc = (0, tt) are the same, so we have 
only three conditions.) 

{At + fif + Al = ifiBH.f, + Al = {piBH,)\ {At - + Al = ifiBH^f. (8) 

When one of these equations ([5]) is satisfied, the energy gap closes. From these conditions, we find that there are 
at least 7 regions of parameter space, which may be topologically distinct, as shown in Fig [TJ We will explore the 
topological numbers associated with these different regions and classify the topological phases of our system in SecfVl 

B. d + id-wave Rashba superconductor 

1. Case of A{k) = A^^^cos ky — cos kx) + i A^f-' sin kx sin ky 

For this d + id-wave gap, the second equation in ^ is met either when Co{k) = or when A(fc) ~ 0. This condition 
is satisfied at fc = (0,0), (7r,0), (0,7r), {tt,tt). Substituting these fcs' in the first equation in ([T]), we have 

(4^ + A^)2^(MBff.)^ ^l'+4{A^J^f^{^lBH,)', {At- ^lf ^ {ii^H^f. (9) 

The conditions obtained here are very similar to those for s-wave pairing. However, there is an important difference 
between them. The first and the last equations of the gap closing condition © do not depend on the amplitude of 
the pairing gap. Because of this feature, even for a relatively weak Zeeman field where the orbital depairing effect is 
negligible, the gap can close, and the topological phase transition occurs, provided that the chemical potential ^ is 
properly tuned as /i ~ ±4i. As will be seen later, this point is crucially important for the feasibility of the realization 
of the non-Abelian topological order in the d -\- id-wave pairing case, compared to the s-wave pairing state. 

2. Case of A{k) = A^^'(sin^ k^ — sin^ ky) + iA^^^ sin k^ sin ky 

As in the previous case, the second equation in ([7]) is met either when Co{k) = or when A(fc) = 0. Accidently, 
both of them are satisfied at the same momenta k = (0, 0), (0, tt), (tt, 0), (tt, tt). Substituting those values into the first 
equation in ([7]), we have 

{4t + ^i)^ ^{fiBH^y, ^^^{pBH^f, {At-fj.)^ ^{(IbH,)^. (10) 

In this case, all of the gap closing conditions (fTO|) do not depend on the pairing gap, and thus, the topological phase 
transition can occur even for a weak magnetic field for /i ^ ±At, 0. The difference between the second equation of © 
and that of ([TO)) yields a slight difference of the parameter regions where a topological order occurs. 
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IV. DUALITY RELATION IN BDG HAMILTONIAN 

As discussed in ref.— , an underlying mechanism of the reahzation of the non-Abehan topological order in the 
Rashba s-wave superconductor is understood in terms of the duality relation satisfied by the model ([T]); i.e. the BdG 
Hamiltonian 'H(fc) is unitary equivalent to the following dual Hamiltonian 'Hj^{k) 

Hu{k) = Dn{k)D^ = ( , f ^(t\r "f"'? "'^'^^^ " r'aS ^ ■ ""0 ' (11) 

^ ' ^ ' yi[e(k) + ilmA{k)\ay + iaCo{k)ay ■ a —ReA(k) + finHzaz J 

where D is the constant unitary matrix given by 

As is shown in SecEl the dual transformation accomplished by the constant unitary matrix does not change the first 
Chern number of the system. Therefore, the original Hamiltonian has the same topological properties as the dual 
one. 

It should be remarked here that the Rashba spin orbit interaction aCo{k) ■ a in the original BdG Hamiltonian 'H(fc) 
induces "the p-wave gap function" — a>Co(fc) ■ o'Cy in the dual BdG Hamiltonian. However, this does not necessary 
means that the topological properties of our system is the same as those of a p-wave superconductor, since Huik) 
has a non standard kinetic term given by ReAs(fc). Nevertheless, we will show in the following sections that the 
topological order similar to a chiral p + ip-wave superconductor emerges under a large Zeeman field. Furthermore, 
the topological order in our system is much more robust than that of a chiral p + ip-wave superconductor. 

In the most part of this paper, we do not use the dual Hamiltonian ■HdI^), but, instead, analyze the original 
Hamiltonian directly. Our analysis using 'H(fc) in this paper confirm the correctness of the argument based on 

the dual Hamiltonian developed in refi^^. 

V. TOPOLOGICAL NUMBERS 

As a conventional long-range order such as magnetic order is characterized by the existence of a nonzero local order 
parameter, a topological phase is also specified by a characteristic quantity similar to an order parameter; this is a 
topological number. In this section, we evaluate the topological number characterizing the non-Abelian topological 
order realized in Rashba spin-singlet superconductors. 

For two-dimensional time-reversal symmetry broken (TRB) superconductors, on which our discussion is focused, an 
important topological number is the TKNN number (equivalent to the first Chern number) /tknn, which is defined 
as follows. Let us consider the BdG equation 

n{k)\cj)„{k)) ^E„{k)\cj)„{k)). (13) 
By using the normalized occupied states, "the gauge field" A[^\k) is defined as 

At\k) = I {Mk)\dkAn{k)). (14) 

Then the TKNN number is given by 

/tknn = -iT- / dh^dkyF^^\k), (15) 

In Jrp2 

where is the first Brillouin zone in the momentum space, J^^^-* (fc) is the "field strength of the gauge field" A[ (fc), 
that is J"(-)(fc) = e'^dk,Aj{k). 

The nonzero TKNN number implies the existence of topological order in the system under consideration. In 
general, the nonzero TKNN number allows both the Abelian topological order, for which there are no non-Abelian 
anyons, and the non-Abelian topological order, which is characterized by the non-Abelian statistics. For the Rashba 
superconductor, the Hamiltonian of which Q is a 4 x 4 matrix, we can calculate /tknn directly from the above 
equations. However, here, we exploit a different method for the evaluation of the topological number, which is 
practically easier to be carried out. Furthermore, this method is quite useful for the elucidation of the realization of 
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the non-Abelian topological order. A key idea of our method is to utilize another topological number specific to the 
Rashba superconductors, which is called the winding number f2i 

The winding number is introduced as follows.— Let us consider the particle-hole symmetry of the BdG Hamiltonian, 

r-H(fc)rt = -n*{-k), (i6) 

where T is given by 

For ky = or tt, it is found that the BdG Hamiltonian ^ of our model satisfies H*{—k) = T-L{k). Thus the particle-hole 
symmetry yields that 

[T,'H{k)]+ = Q (18) 
for ky = 0, TT. From this relation, it is found that if we take the basis where F has the diagonal form as 

r=(' Y -t ^"l' (19) 



2x2 



then T-L{k) at ky tt becomes off-diagonal, 



q(k) 



By using qik) in the above, the topological number I{ky) is defined as 

Hky) = ^ y dk^ir [q~Hk)dkMk) - q^-\k)dkj{k)] . (21) 

We call I{ky) as the winding number in the following. 

As is shown in Appendix IB 31 the winding number and the TKNN number satisfies 

(-_2^^/tknn ^ ^_^y{o)-i{7T) ^ ^22) 

Therefore, from the winding number, we can determine (— 1)^tknn^ ^j^g index (— i)^tknn jg of particular interest since 
it give a hallmark of the non-Abelian topological phase: when (— 1)^tknn = there are an odd number of Majorana 
zero modes in a vortex, which implies the vortex is a non-Abelian anyon. In the following subsections, we calculate 
the winding number and the TKNN number for the cases of the s-wave pairing state and the d + id pairing state. 
The obtained phase diagrams for the non-Abelian topological phases are summarized in Fig[TJ 

A. s-wave pairing 

We consider an s-wave NCS superconductor with the gap function 

A(fc) = A,. (23) 
For the s-wave NCS superconductor, q{k) is given by 

lik) = - [£{k) - Mb^^zCt^ - asiiikxay] + lAgay, (24) 

thus its determinant is 

detq{k) = e{k)'^ ~ [iiYiUzf - o? sin^ Ai^j; + - 2iaAs sin A:^;. (25) 
Denoting the real (imaginary) part of det(7(A:) as toi(A:) (mi(ki)\ we have 

TOi(fc) = e(fc)2 - (^B-ffz)^ - a^sin^ A:^; + A^, m2(fc) = -2aAs sin Ai^,. (26) 



7 






FIG. 1: The diagrams of the non-Abehan topological phase for spin-singlet NCS superconductors, (a) s-wave case, (b) 
ci + id-wave case with A(fc) — A^^\coaky —coakx) + i A^^^ sin kx sin ky. (c) d 4- id-wave case with A(fc) = A^^' (sin^ fcj, — 
sin"^ ky) + iA'^^'' sin kx sin ky . The topological numbers ((— l)"'^''; 7(0), /(tt)) are given only in the phases supporting a non- 
Abelian topological order. In each case, there are four different non-Abelian topological phases. In the s-wave case, the 
non-Abelian topological phase is realized only when the Zeeman magnetic field satisfies (fj-sHz)^ > (As)^, but in the d -I- id- 
wave cases, the non-Abelian topological phase can be realized even for a small but nonzero H^- 



From the formula (|B22[) in Appendix IB 21 the winding number is evaluated as 

Hky) - I [-sgn [£(0, ky f - {fiBH^ f + Aj] + sgn [e(^, fc,)^ - {^^^H^ f + A^] ] . (27) 

In Table m we summarize the winding number I{ky) calculated from this equation. We also list ( — 1)^tknn obtained 
from the formula (P^ . Note that the topological numbers can change only when one of the gap closing conditions (|S]) 
is met. From Table HI we obtain the diagram of the non-Abelian topological phase for s-wave pairing shown in Fig[T] 
(a). 

To get a better understanding of the origin of the non-Abelian topological order, we, here, present some physical 
discussions about the phase diagram shown in FigH] (a). We consider two different but complementary arguments. 
The first one is based on mapping from the s-wave pairing state to an effective spinless p-wave pairing state in the 
chirality basis. The second one is the argument based on the duality relation introduced in Sec lIVI The former is 
applicable to the case of HbHz ^ A^, while the latter is particularly useful in the vicinity of the topological phase 
transition point hbHz ^ As. In this sense, these two arguments are complementary. 

We, first, present the first argument which utilizes the mapping onto an effective p-wave pairing state in the chirality 
basis representation, and is applicable to the case of /ie^z 3> A^. In the chirality basis, the SO coupling term and 
the Zeeman coupling one in the Hamiltonian are diagonalized, leading to the two SO split bands. As is shown in 
Appendix [Cl in the parameter region where the non-Abelian topological phase is realized, we can map our model 
dH) into a spinless chiral p-wave superconductor in the chirality basis, provided that HbHz ^ A^. For /i < 0, the 
low-energy effective Hamiltonian in this case is, 

~ / e{k)~/Ae{k) {a£oAk)-iaCov{k)){^{k)/Ae{k)) 

"'-^''> - [ (aCoAk) + iaCay{k)){A*{k)/Ae{k)) -e(fc) -I- Ae(fc) 
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and for /i > 0, it is 

^ s{k) + Ae(k) (aCoAk) + iaL^y{k)){A{k) / As{k)) 

"■+^^> " \ {aCo^k) - taCoyik)){A*{k)/Ae{k)) -e(fc) - Ae{k) 

with Ae(fc) — {aCo{k))^ + (/ie-ffz)^- For the s-wave pairing state in the original Hamihonian, the gap function in 
(f28| or ^ is given by 

{aCox{k)T'ict£oy{k)){A{k)/Ae{k)) ia{sinkx ^ ismky){As/ fiBHz), (30) 

thus, for both /i's, the chiral p + ip wave superconductors are reahzed in the chirahty basis. By using the effective 
Hamihonian, the phase diagram of our system can be understood more intuitively. The above effective Hamiltonian 
([28| or ([29l) is obtained by using the fact that, when hbHz ^ Ag, and fi is in the region of the phase diagram where 
the non-Abelian topological order is realized, only one of the two Fermi surfaces survives and the other is pushed 
away by the Zeeman magnetic field. As a result, the spinless chiral p + zp-wave superconductor is realized effectively. 
The non-Abelian topological phase obtained here is effectively the same as that in the spinless chiral p + zp-wave 
superconductor. 

We, now, present the second argument on the origin of the topological order which is based on the duality relation)^ 
To grasp the physics shown in the phase diagram, let us see what happens at the transition between the trivial phase 
{i.e. the phase with hbHz — 0) and the non-Abelian topological one ((— l)-fTKNN _ —l^. From Fig[T] (a), it is found 
that if we increase the Zeeman magnetic field, such a phase transition occurs at fi — ±4t when /ib^^z = As. For 
simplicity, we consider the case with /i = — 4t, where the Fermi surface is close to the F point. A similar analysis is 
possible for the transition at /i = At. In the former case, the gap of the system closes a.t k ~ (0, 0). 

To examine the topological phase transition, it is convenient to use the dual Hamiltonian instead of the original 
one. Using the duality transformation (jlip . we recast the original BdG Hamiltonian into its unitary equivalent dual 
Hamiltonian 'HD[k). Then, we find that for /i ^ — 4i, the dual Hamiltonian 'Hj){k) around k = (0,0) is decomposed 
into the following two 2x2 matrices, 

- ( V ^""-^^ '''i ) ' (31) 

where f and I denote the spin in the basis of the dual Hamiltonian. 

We notice here that these Hamiltonians have a close similarity to the Hamiltonian of the spinless chiral p + ip 
superconductor discussed in^. The spinless chiral p + ip-wave superconductor shows a phase transition between the 
non-Abelian topological phase (or weak-pairing phase in the terminology used in^) and the topologically trivial phase 
(strong-pairing phase) , and the transition is described by the following low energy effective Hamiltonian 

"W (h\ — ( ^p{kx + iky) \ ,„„s 

Hp+^p(k) - ^ _ ^^^^ J , [66) 

where /ip and Ap are the chemical potential and the paring amplitude for the p -I- ip-wave superconductor, respectively. 
For Hp > 0, the state is topologically trivial, and for /ip < 0, the state supports non-Abelian topological order. The 
phase transition occurs at /ip = 0. By identifying A^ — ^bHz and a in with //p and Ap in (I33p . respectively, the 
similarity between 'H-^-[[k) and Jip+ip^k) is evident. This similarity immediately implies that the phase transition at 
HbHz — As is also accompanied with the emergence of the non-Abelian topological order. (We also have a similar 
phase transition at ^ibHz = —As by decreasing the Zeeman magnetic field. From the similarity between H^^k) and 
T~ip+ip{k), this transition is also found to be accompanied with the emergence of the non-Abelian topological order.) 

Indeed, using the dual Hamiltonian (PT|) . we can see directly that the TKNN number change by A/tknn = — 1 at 
the transition. For this purpose, we slightly generalize the gauge field A| ''(fc) as follows, 

a{-)(u\ _ \ 'i'T.E„<o'^^n{k)\dk^(t)n{k)) for/i = x,?/, 

^ ' ~ I *Ei;„<o('^"(fc)|5/.B//>n(fc)) for /^ = ^^^^ 

where the z-component is introduced additionally. Then, consider the rectangle illustrated in Figl2l Here the top 
face BZ(I) and the bottom one BZ(II) denote the first Brillouin zones of the system after [hbUz > As) and before 
[ixbHz < As) the transition, respectively. Because of the periodicity of the Block wave function (/)„(fc) in the 
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momentum space, "the magnetic field" \k) on the side faces Si (i = 1,2) are identical to that on the opposite 
ones S[. Therefore, the change of the TKNN number 



TKNN 



^[ dk^dkyT^-\k) - ^ [ dk.,dkyT^-\k). 

JBZ(I) JBZ{II) 



(35) 



5(1) ^" ^BZ(II) 

is rewritten as the total "magnetic field" penetrating the surface of the rectangle dV , 



A/tknn = 77-/ dk^dkyT'^-\k) - ^ [ dk^dkyF'^-\k) 

^''^ JBZ(I) JBZ(II) 

+ — [ dkydk:,F^-\k) - — I dkydk,F'^-\k) 
27r Js, 27r Jg, 

+ _L /■ dk,dk,T''-\k) - ^ [ dk,dk,T^-\k) 

= J_ /■ d5J-(-)(fc), 
27r Jay 



(36) 



with /c^ = hbHz — Ag. Therefore, A/tknn is nonzero only if the "magnetic monopole" of \k) exists inside the 
rectangle: If there is no such a source of the magnetic field, A/tknn should be zero from the Gauss-Bonnet theorem. 
On the other hand, if the magnetic monopole exists, we have a net magnetic fiux penetrating the surface of the 
rectangle. 







BZ(I) 










Si 


BZ(n) 


/ 



FIG. 2: Topological phase transition at hbHz — As. BZ(I) and BZ(II) indicate the Brillouin zones after and the before the 
transition. The gap closes at (k^, ky) = (0, 0) when fisHz = As. S,; and (i — 1, 2) are the side faces of the rectangle in which 
the top and bottom faces are BZ(I) and BZ(II). For simplicity, only Si and S'l are explicity indicated. 



Actually, we have a magnetic monopole located at {kx, ky,kz{= /xei/z — A^)) = (0, 0, 0), where the gap of the system 
closes. The magnetic charge can be read from the dual Hamiltonian H^^{k). By rewriting H^^ik) as H^^{k) = R{k) ■ a 
with {Ri{k) , R2{k) , R^ik)) = (aky,—akx,—kz), the monopole charge is given by, 

Q ^ ^ I dSe^,R-{dk^^Rxdk^R) (37) 

where R{k) = R{k)/\R{k)\ and 5^ a small sphere surrounding the gap-closing point. Noting that the right hand side 
of (j37p counts the number of times the unit vector R wraps around the origin, we obtain Q — —1. Therefore, from 
P6p . we immediately find that A/tknn = — 1- Before the transition, the system is topologically equivalent to the 
ordinary s-wave superconductor without the Zeeman magnetic field, thus /tknn = 0. Therefore, we have /tknn = — 1 
after the transition)^ Again, this result indicates that the system after the transition belongs to the same topological 
class as the spinless chiral p + ip-wave superconductor with /tknn = — 1- 

As mentioned before, the above two arguments are applicable, respectively, to the different parameter regions, and 
thus, they are complementary. It is noted that these arguments are also straightforwardly applied to the case of the 
d + id pairing state discussed in the next subsections. 

As seen from Table HI the non-Abelian topological order (i.e. (— 1)^tknn — _2) appears only in the case that the 
Zeeman energy iibHz is larger than the superconducting gap A.,. As is well-known, the Rashba superconductors are 
stable against the Pauli depairing effect due to applied magnetic fields even for /iB^^z > A^ when the magnetic field is 
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a) M < -2t 





(- 




1(0) 


JW 






1 








{At + + < {fijiH. f < + A? 




-1 


1 





/i^ + A2 < (^iB//,)2 < (4t - + A'i 




-1 





1 


(4t - + A? < (mbH.)' 




1 








b) < < 




(- 




1(0) 


7W 


0<(AtBJYz)'<M' + A^ 




1 








+ A^ < ifiBH^f < {4t + + Al 




1 


-1 


1 


(At + iif + Al< i/iBH^f < {4t - + Al 




-1 





1 


(4t-M)2 + A? < {psH^ f 




1 








c) < fi<2t 




(- 




1(0) 


/W 


< ifiBH. f < fi^ + Al 




1 








^? + Al< i/iBH^ f < {4t - tif + A^ 




1 


-1 


1 


{4t - + Al< {pBH. f < {4t + fif + Al 




-1 


-1 





{4t + +A'i< (mbH,)2 




1 








d) 2t < n 


{flBH.f 


(- 




1(0) 


I(7V) 


Q<{^JiBH..f < {4t~^if + Al 




1 








{4t - fif + A? < iiiBH. f < At' + A? 




-1 





-1 


+ A? < {fiBH.f < {4t + fif + A'i 




-1 


-1 





{4t + fif + A? < ifiBH.f 




1 









TABLE I: The TKNN integer /tknn and the winding number I(ky) for 2D s-wave superconductors with the Rashba couphng. 
(_1)^TKNN _ corresponds to the non-Abehan topological phase. 

applied perpendicular to the xy-plane^^^i^ as long as the Rashba SO interaction is sufBciently strong. However, there 
is also the orbital depairing effect due to applied magnetic field. An important question is how the superconductivity 
survives the orbital depairing effect for such a strong magnetic field pbHz > A^. One possible scenario is to realize this 
system in the proximity between a superconductor and a semiconductor as proposed in refs.^^"^^ Another possibility 
is to realize it in strongly correlated electron systems for which the orbital depairing field is large. Also, one more 
promising scheme is to utilize ultracold fermionic atom as proposed in2&. This issue will be discussed in more detail 
in Secim 

B. d + id wave pairing 

1. Case of A(k) — A^^^ (cos ky — cos kx) + iA^P sin k^ sin ky 
For the case of the d + zci-wave superconductor with the gap function 

A(fc) = A^^-* (cos ky — cos kx) + i^^^'^ sin sin ky, (38) 
q[k) and its determinant are given by 

q(k) = — [s(k) — ^j,bHz<Jz — asink^ay] + iA'"^\cosky — cos kx)ay, (39) 
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detq{k) = — (ij-bHz)^ — o? sin^ + ( A^^^'' (cos fcj^ — cosfc^;))^ — ^ia'smk^l^^^ {co%ky — cosfcx)- (40) 

Therefore, the real and imaginary parts of the determinants are 

■mi{k) — e{kY — {(J-eHz)^ ~ c? sin^ k^ + (A^"*^-* (cos fcj^ — cosfc^))^, rai(k) = —2iasinkx/S.^p [cosky — coskx). (41) 

To apply the formula (|B22[) . we slightly change m2{k) as — 2iasinfc2:A^^''(cosfcy — cosfc^;) + 6 {6 ^ 1). After using 
the formula (jB22l) . we put (5 = again. Since I{ky) is a topological number, this procedure does not change the value 
of I{ky). In this manner, we estimate the winding number I{ky) as 

^(0) = I [-sgn[£(0, 0)2 - ifXBHzf] + sgn[£(^, 0)^ - (/iBi/.)^ + 4(A|,^))2]] , (42) 

and 

H^) = I [-sgn[£(^, 7r)2 - itiBHzf] + sgn[e(0, tt)^ - [tiBH^f + 4(A(^V]] • (43) 

We summarize the winding number I{ky) and (— i)^tknn jj^ Table HIl The non-Abelian phases are realized in the 
parameter regions that (— 1)^tkkn _ _i pYom Table II, we obtain the phase diagram shown in FigH] (b). 

In a manner similar to the NCS s-wave superconductor, if ^bHz 3> A(fc), the obtained phase diagram can be 
understood by the effective Hamiltonian p8| for /i < or ((29|) for /i > obtained in the chirality basis. For 
simplicity, suppose that A^^-* = l^'^p = A^. The gap function in ((^5)) or (|29l) yields 

{aCQx{k) ± ia£Qy{k j){A{k) / Ae{k)) ^ ^ia{smkx ± i smky){cosky — cosk^ + i sinfej, sin ky) {A ^ / fiB Hz), (44) 

thus the chiral / + z/-wave superconductor or the chiral p + ip-wave superconductor is realized both in the effective 
Hamiltonians T-L^{k). (Whether the p-wave state or the /-wave state realizes depends on the relative chirality between 
the d + id order parameter and the p-wave factor in (|44p which stems from the SO interaction.) Therefore, the non- 
Abelian topological phases obtained here are effectively the same as that in either the spinless chiral / + i/-wave 
superconductor or the spinless chiral p + ip-wave superconductor. 

One remarkable point observed from Table HIl (or Figlljb)) is that, in contrast to the s-wave pairing case, for 
the d + id-wave pairing state, the non-Abelian topological order appears even for small but nonzero magnetic fields, 
provided that fi ~ ±4t. Thus, in this case, we do not need to worry about the orbital depairing effect due to 
applied magnetic fields. This point makes it easier to realize the Majorana fermion state in the d + id Rashba 
superconductor than in the s-wave pairing state from the perspective of the stability against applied magnetic fields, 
though, unfortunately, the experimental realization of d 4- id superconductors has not yet been established to this 
date. 

Since the d + id pairing state breaks time-reversal symmetry, the TKNN number is nonzero even for zero magnetic 
fields. In the absence of the Zeeman field, we can evaluate the TKNN number /tknn directly. In this case, we 
smoothly eliminate the Rashba SO interaction by setting a without gap closing. This means that the TKNN 
number of the d -\- id-wave NCS superconductor is the same as that of the d-\- id superconductor without the Rashba 
coupling. 

2. Case of A(fc) = A^^^ (sin^ kx — sin^ ky) + i^^P sin kx sin ky 

For the d -\- id-wave NCS superconductor with the gap function 

A(fc) = A^^^ (sin^ kx - sin^ ky) + iA^^ sin k^ sin ky, (45) 

q{k) is given by 

q{k) = - [e{k) - ^bHzCTz - asinkxcry] + iA^'p sin^ k^ay. (46) 
Thus its determinant becomes 

detq(fc) = e(fc)2 - (i^bHz)^ - sin^ k^ + (A^)^ sin"* k^ - 2iaAP sin^ k^, (47) 
and the real and imaginary parts of the determinant are 

TOi(fc) =e(fc)2 - (hbHz)^ - si-a^ kx + iAP)^sm^kx, m2{k) = -2aA|;^^ sin^ fc^. (48) 
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a) fi< -2t + {A['-^) 


y2t 










1(0) 




0< (mbH.)' < (4t + M)' 


1 








(4t + M)'<(MBi/.)'<M'+4(Aii')' 


-1 


1 





^? + 4(A(i')2 < (/iBi/.)' < (4t - ^lf 


-1 





-1 




1 








b) -2t + {A^^^f/2t < fi<0 






7(0) 




0<(MBi/.)^<M'+4(AW)^ 


1 








/i' + 4(A('^)2 < (mb//.)' < (4t + m)' 


1 


-1 


-1 


(4f + fif < {/iBH^f < {4t - t^f 


-1 





-1 


(4t-/i)2 < ifisH^f 


1 








c) 0< fi<2t- (AJ;^' 


r/2t 






(mbH.)' 




1(0) 




0< (mb^.)' <m'+4(AW)' 


1 








/i' + 4(A(i))2 < (;.b//.)' < (4t - m)' 


1 


-1 


-1 


(4f - fif < {iiBH.f < {At + ^lf 


-1 


-1 







1 








d) (A^'')V2t < 


(mbH.)' 




1(0) 


/W 


< (MBi?.)" < (4t - m)' 


1 








(4t - fif < ifiBH.f < + 4(A<^')2 


-1 





1 


+ 4(A(i')2 < (^ei/.)" < (4t + m)' 


-1 


-1 





(4t + M)' < ifiBH,)" 


1 









TABLE II: The TKNN integer /tknn and the winding number I{ky) for the 2D d + id-wave NCS superconductor with 
A{k) — A^^' (cos ky — cos kx) + iA^^' sin kx sin ky. (— 1)^tknn — _i corresponds to the non-Abelian topological phase. 

In a manner similar to the previous d-wave case, we regulate ni2{k) as ni2 (k) -2aA^^^ sin^ + (5 < 1). Then, 
we obtain 

HK) = \ [-sgn [£(0> hf - iP^BH.f] + sgn [e(^, kyf - (Mei?.)']] • (49) 

We summarize the winding number I{ky) and (— 1)^tknn jj^ Table Hill The results are also summarized in the phase 
diagram FigU] (c). 

When fiBHz ^ ^{k), the obtained non-Abelian topological phase can be understood by using the effective Hamilto- 
nian ([25)1 or (|29l) in a manner similar to the previous cases. In the parameter region where the non-Abelian topological 
phase is realized, only one of the Fermi surfaces survives and the other is pushed away by the Zeeman magnetic field, 
then, it is found that the non-Abelian topological phases obtained here are effectively the same as that in cither the 
spinless chiral / -I- if -wave superconductor or the spinless chiral p + ip-wave superconductor. 

As in the previous subsection, the non-Abelian topological order is realized even for small but nonzero magnetic 
fields, provided that fi ^ ±41. This is a generic feature of the d -\- id-wave pairing state. However, in contrast to 
the case of A(fc) — A^^-* (cos ky — cos fc^) + i^li^ sin k^ sin ky , the parameters which distinguish different phases do not 
depend on the superconducting gap function. 
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a) M < -2t 









7(0) 


7(7r) 


< {^^BH,f < {4t + fif 




1 








{4t + fif < i^iBH.f < 




-1 


1 





< {f^BH. f < {4t - 




-1 





1 


(4t-/i)2 < ifisH.f 




1 








h) -2t< n<0 




(- 




7(0) 


7(7r) 


< ifiBH.f < 




1 








< il^BH^f < {4t + tif 




1 


-1 


1 


{4t + fif < {/iBH^f < {4t - t^f 




-1 





1 


i4t-nf < {fiBH.f 




1 








c) < fi<2t 


{flBH.f 


(- 




7(0) 


lin) 


< if^BH^f < ^2 




1 








< if^BH^ f < {4t - tif 




1 


-1 


1 


{4t - fif < (nBH. f < {4t + fif 




-1 


-1 





{4t + ^f < ifiBH.f 




1 








d) 2t < n 








7(0) 




< {^^BH,f < {4t - m)' 




1 








{4t-fif < iiiBH.f < ti^ 




-1 





-1 


< i^^BH,f < {4t + fif 




-1 


-1 





{4t + iif < ifiBH,)" 




1 









TABLE III: The TKNN integer 7tknn and the winding number I{ky) for the 2D d + id-wave NCS superconductor with 
A(fc) = AJ^^' (sin^ kx — sin^ ky) + i^f ^ sin k^ sin ky. ( — 1)^tknn _ corresponds to the non-Abelian topological phase. 



VI. MAJORANA CHIRAL EDGE STATE 



In this section, we investigate edge states for the 2D spin-singlet NCS superconductors numerically. From the 
bulk-edge correspondence, a non-trivial bulk topological number implies the existence of gapless edge states. In the 
case of TRB superconductors, the gapless edge states is a chiral Majorana fermion mode. We confirm this in the 
following. 



A. s-wave 



To study edge states, we consider the lattice version of Hamiltonian ([T]). For an s-wave NCS superconductor, the 
lattice Hamiltonian is given by 



(50) 
(51) 

(52) 
(53) 
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where i = {ix,iy) denotes a site on the square lattice, c]^ {cia) the creation (annihilation) operator of an electron 
with spin a at site i, and A = a/2. The sum j) taken between the nearest-neighbor sites. Suppose that the 
system has two open boundary edges ed ix — and ix — Nx, and impose the periodic boundary condition in the 
y-direction. By solving numerically the energy spectrum as a function of the momentum ky in the j/-direction, we 
study edge states. 

We illustrate the energy spectra for the 2D s-wave NCS superconductor with edges at ix = and ix — 30 in Figs 
[3]and[ll In Fig. [3] (Fig. H]), ^ satisfies < —2t {—2t < fi < 0), and the corresponding bulk topological numbers are 
given in Table |T] (a) (Table |T] (b)). We find that if the TKNN number is odd, odd numbers of gapless edge modes 
appear. In this case, the non-Abelian topological order appears, and the edge zero mode is a chiral Majorana fermion. 
It is also found that when the winding number I{ky) is non-zero for ky = 0, or tt, the energy of the gapless edge mode 
becomes zero at this value of ky. Therefore, all the results are consistent with the existence of the correspondence 
between the bulk topological numbers and the edge spectra. 




Ill E l\ E 




FIG. 3; The energy spectra of the 2D s-wave NCS superconductor with open edges at i^: = and ix — 30 for n < —2t. Here 
ky denotes the momentum in the y-direction, and ky £ [— tt, tt]. We take t = 1, /i = —2.5, A = 0.5, and As = 1. The Zeeman 
magnetic field is (I) fisH, = 0, (II) fisH, = 2, (III) hbH, = 3, and (IV) pB-ffz = 7. The cases of (I), (II), (III), and (IV) 
correspond to, respectively, the four regions in Table Ha). The non-Abelian phases are (II) and (III). The two gapless modes 
found in (II) correspond to edge modes for two open edges, respectively. Thus, they are chiral. The same is also true for the 
gapless modes in (III). 



B. d + id wave 



1. Case of A{k) = A^J'\cos ky — cos k^) + iA''^' sin k^ sink 



(2), 



The lattice Hamiltonian for the 2D d + id-wave NCS superconductor with A(fc) — A^^''(cosfcy — cosfea,) 

(2) 

iA^ sin kx sin ky is given by "H = 'Hkin + "Hso + "Hs with 



(1) 



A 



(2) 



'4+x^'^i-l + ^l-x^'^i-l ^l+yt'4l '^i-yt'^l-l 



"i+x+yt'^'il ~^ '^i-x-y-\'^ii '^i+x-y^^ii 



H.c. 



(54) 
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FIG. 4: The energy spectra of the 2D s-wave NCS superconductor with open edges at ia; = and ~ 30 for —2t < fj. < 0. 
Here ky denotes the momentum in the y-direction, and ky G [— tt, tt]. We take t = 1, fi = —1, X = 0.5, and As — 1. The Zeeman 
magnetic field is (I) fiaH, = 0, (II) fisH, = 2, (III) ^bH, = 4, and (IV) pe-ffz = 6. The cases of (I), (II), (III), and (IV) 
correspond to, respectively, the four regions in Table Ulb). The non-Abelian phase is (III). 



The kinetic term "Hkin and the Rashba SO interaction Hso are the same as (|5ip and ([52|). respectively. 

In a manner similar to the s-wave NCS superconductor, we obtain the energy spectra for the system with edges 
at ij; = and ix = 30 numerically. We illustrate the energy spectra for the 2D d-wave NCS superconductor with 



the gap function A(fc) 



^d\cOa ky 



coskx) + jAJj sinfca; sinfcj^ in Figs [5] and [6] In Fig. [5] (Fig. ^ satisfies 



H < 



-2t+ {ls!]pY /2t (— 2i-|- {ls!]pY /2t < fi < 0), and the corresponding bulk topological numbers are given in Table 



im (a) (Table In] (b)). We find that if the TKNN number is odd, odd numbers of gapless edge modes appear. The 
non-Abelian topological order (and hence the chiral Majorana fermion mode) is realized in this case. In addition, it 
is found that if the winding number I{ky) [ky = 0,7r) is non-zero for some ky, the energy of the gapless edge state 
becomes zero at this value of ky. These results are consistent with the bulk-edge correspondence. 

Since the d + id-wave superconductor breaks time-reversal symmetry even in the absence of a magnetic field, the 
TKNN number is nonzero for = 0, and there are four chiral edge modes, which are seen in FigsE] (I), E] (I), [7] (I), 
and [5] (I). Since there are even numbers of gapless edge modes, they do not behave as non-Abelian anyons. These 
edge modes crosses the zero energy at some ky ^ 0,±7r. In contrast, the Majorana fermion mode associated with 
the non-Abelian topological order crosses the zero energy at ky = or tt corresponding to the nonzero values of the 
winding number I{ky) at these points. 



2. Case of A(fc) = (sin^ k^ — sin^ ky) + iAjj sin k^ sin ky 
In this case, we use the lattice Hamiltonian % = "Hkin + ^so + with 



A 



(1) 



.(2) 



4+2£t'^4 + '^i-2x\'^ii "^i+2yt'^4 '^i-2j/t''4 



'^i+x+y-\'^ii + '^i-x-y^'^ii '^i+x-y^'^ii ^i-x+yt'^il 



H.c. 



(55) 



Here, the kinetic term Hi^in and the Rashba SO interaction Hso are the same as and ([5^ . respectively. 



16 



1 £ n £ 




111 E IV £ 




FIG. 5: The energy spectra of the 2D d + id-w&ve NCS superconductor (A(fc) — AJj^' (cos ky — cos kx) + i^^P sin fci, sin ky) with 
open edges at i^, = and ix ~ 30 for ^ < — 2i; + A^^-'^/2t. Here fcj, denotes the momentum in the y-direction, and ky £ [— vr, tt]. 
We take t=l, fj,= -2.5, A = 0.5, A^;^^ = 0.5, and Af'' = 0.8. The Zeeman magnetic field is (I) fisH^ = 0, (II) fisH^ = 2, 
(III) fiBHz = 3, and (IV) I-IbHz = 7. The cases of (I), (II), (III), and (IV) correspond to, respectively, the four regions in Table 
|TTla). The non-Abelian phases are (II) and (III). 



We calculate the energy spectra for the system with edges at z^; = and ix — 30 numerically. The energy spectra 
for the 2D d-wave NCS superconductor with the gap function A{k) — A|^^''(sin^ kx — sin^ ky) + iA^ sinfca; sinfcy are 
shown in Figs [7] and [S] In Fig. [7] (Fig. [5]), fi satisfies /i < —2t {—2t < ^ < 0), and the corresponding bulk topological 
numbers are given in Table lllll (a) (Table Hill (b)). When the TKNN number is odd, odd numbers of gapless edge states 
appear, signifying the non-Abelian topological order. It is also found that if the winding number I{ky) {ky = 0,7r) 
is non-zero for some ky, the energy of the gapless edge state becomes zero at this value of ky. These results are 
consistent with the bulk-edge correspondence again. 

VII. MAJORANA ZERO MODE IN A VORTEX OF AN s-WAVE RASHBA SUPERCONDUCTOR 

In this section, we discuss a Majorana fermion mode in a vortex core of an s-wave Rashba superconductor, which is 
one of the important features of the non-Abelian topological order, and relevant to the application to the topological 
quantum computation. We also discuss that, in certain parameter region, the Majorana mode in the s-wave Rashba 
superconductors is strongly stable against thermal noise and inter-vortex tunneling, which are serious obstructions to 
the implementation of the topological quantum computation in the case of p + ip superconductors. 

A. Majorana solution of the Bogoliubov-de Gennes equation 

The non-Abelian topological order is characterized by the existence of the Majorana zero energy mode in a vortex 
core, which leads to the realization of the non-Abelian statistics. In this section, we demonstrate that the Majorana 
fermion mode in a vortex core exists for an s-wave Rashba superconductor when the Zeeman energy (IbHz is larger 
than the superconducting gap A, by solving the Bogoliubov-de Gennes (BdG) equation for a single vortex of the 
superconducting order parameter; A(r) = Aexp{in9) with n vorticity. In ref.^^, when n is odd, the Majorana bound 
state for a vortex of the SO interaction is obtained by the analysis of the BdG equation for the dual Hamiltonian 
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FIG. 6: The energy spectra of the 2D d + id-wave NCS superconductor (A(fc) = AJj (cosfcy — cosfc^;) + iA]j sinfca; sinfcy) 
with open edges at ix — and ix = 30 for —2t + A^^^^ /2t < fi < 0. Here ky denotes the momentum in the y-direction, and 
ky e [-n,TT]. We take t = l, = -1, X = 0.5, A^^' = 0.5, and A^^ = 0.8. The Zeeman magnetic field is (I) fisH^ = 0, (II) 
fJ'sHz = 1.6, (III) I-IbHz = 3.1, and (IV) fisHz = 6. The cases of (I), (II), (III), and (IV) correspond to, respectively, the four 
regions in Table Ulll b). The non-Abelian phase is (III). 



PT|) . The existence of the Majorana bound state in a vortex of the SO interaction strongly implies that there exists 
the Majorana bound state also in a vortex of the superconducting condensate, because these two vortex states are 
related by a singular gauge transformation. In the following, we will find the zero energy Majorana solution for a 
vortex core of the superconducting condensate explicitly in the case that the vorticity n is odd, and the condition 
fisHz > A > is satisfied. This result is in accordance with the finding of Sau et al<21 In the following analysis, 
instead of using the dual Hamiltonian (fTTj) . we deal with the BdG equation for the original Hamiltonian (j4]) directly. 
Since we can not obtain the exact solution of the BdG equation analytically, we adopt the following approximation 
scheme. As clarified in the previous sections, the non-Abelian topological order appears when the Fermi level crosses 
A:-points in the vicinity of the F point or the M point in the Brillouin zone. In this situation, there are two Fermi 
surfaces split by the Rashba SO interaction. One is located in the vicinity of the F (or M) point with the Fermi 
momentum kp ^ (or (7r,7r)), and the other has the large Fermi momentum kp ^ 0, (tt, tt). Since the momentum is 
not a good quantum number in the presence of a vortex core, a bound state in the vortex core generally consists of 
a superposition of quasiparticles from both of these two Fermi surfaces. However, according to the discussion based 
on the duality Hamiltonian given in Sec IV A[ it is strongly suggested that in the vicinity of the topological phase 
transition point (J-bHz ~ A, quasiparticles with fc ^ (or (tTjTt)) play a very important role in the realization of the 
non-Abelian topological order. This implies that when the Zeeman energy is close to the superconducting gap, the 
Majorana fermion mode is mainly formed by quasiparticles with A:f ~ or (tt, tt) rather than those with kp ^ 0, (tt, tt), 
in the long-distance asymptotic regime away from the center of the vortex core. Thus, in the following, we try to 
construct an approximate solution for the zero energy mode in a vortex from quasiparticles with A:^^ or (tt, tt). As 
will be shown below, the vortex core for this approximated solution has a characteristic length ~ vp /{(J-eHz — A), 
while quasiparticles with kp 0, (tTjTt) give contributions to the vortex core bound state with a characteristic length 
~ vp/A. Thus, the approximated solution presented in the following is valid when < (J-bHz — A < A in the 
long-distance asymptotic regime. 

To solve the BdG equation, we choose the gauge for which the gap function is real by applying the gauge transfor- 
mation eA eA — nW9/2, A exp{in9) — > A. Then, the BdG equation is 



(56) 
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FIG. 7: The energy spectra of the 2D d + id-wave NCS superconductor (A(fc) = A^^'(sin'^ fc^ — sin"^ ky) + iA^^' smfc^; sinfcy) 
with open edges at ia; = and ix ~ 30 for ^ < —2t. Here ky denotes the momentum in the j/-direction, and ky € [— 7r,7r]. We 
take t = l,ij,= -2.5, A = 0.5, A^^' = 1, and A^ ' = 1. The Zeeman magnetic field is (I) fisH^ = 0, (II) fiaH^ = 2, (III) 
fJ^sHz = 3.5, and (IV) ^lbH^ — 7. The cases of (I), (II), (III), and (IV) correspond to, respectively, the four regions in Table 
lllll a). The non-Abelian phases are (II) and (III). 



" -l\iuy -e{k + eA- 'jVe) + g{k + eA~ ^Ve) -a* + ha, J ' ^ ^ 

Here ^'^ = (u^, u^, w;), e{k) = ^ - /i, g(k) = 2X{ky, -k^,0), k = -iV , and h = ^bH,. Assuming < Hc2, we 

neglect eA compared to n\76/2 in e{k — eA + f V0) and g{k — eA+^\79)M We also assume that A = for r < Tc, and 
A 7^ for r > Tc, where Tc ^ ^. For simplicity, we consider the case of /i = 0, for which one of the two SO split bands 
crosses the T point fc = 0. This is a typical situation which realizes the non-Abelian topological order (and hence 
the Majorana zero mode) as discussed in the previous sections. In the following, we restrict our analysis to the zero 
energy state with E = 0. Furthermore, we impose the condition that the magnitude of the SO interaction is much 
larger than the Zeeman energy scale; i.e. h <gC mX^. Under this condition, we can find the following approximated 
solution for the zero energy mode in a vortex core with the odd vorticity n. For r < rc, 

Mr,e)^A^e-^-^H'^^i^'^l (58) 

u^(r,e)=-^A^e^-^H^,{^'^), (59) 

and for r > Tc, 

u-f{r,e) = V2Ate-^^e^''*'^i/S3(«^^^^r-), (60) 

2 A 

u^{r,0) = -tV2A-fe'h^^'^'''^H':^l{i'^—^r), (61) 

2 A 

where H^^"^ (z) is the first Hankel function, and the constant A^ is determined by the normalization condition. Also, 
Vaif, 9) = Ua-{r, 9). The solution for r > Vc and those for r < Tc can be matched at r = Vc by using the asymptotic 
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FIG. 8: The energy spectra of the 2D d + id-wave NCS superconductor (A(fc) = A^^-'(sin^ ~ sin^ ky) + iA^'^' sinfca; sinfcy) 
with open edges at io; = and ix — 30 for ~2t < fi < 0. Here ky denotes the momentum in the y-direction, and ky G [— 7r,7r]. 
We take t = 1, = -1, X = 0.5, AJ^^' = 1, and A^^' = 1. The Zeeman magnetic field is (I) fisH^ = 0, (II) /XB-ffz = 2, (III) 
fJ,BHz = 3.5, and (IV) ^bHz — 6. The cases of (I), (II), (III), and (IV) correspond to, respectively, the four regions in Table 
Ulllb'). The non-Abelian phase is (III). 



form of the Hankel function, as explained in the Appendix |Dl Since there is only one zero energy mode, the above 
solution indicates that there is a Majorana zero energy mode in a vortex core with odd vorticity. Note that this 
approximated solution is constructed from quasiparticles in the vicinity of the F point fc ~ 0. As seen from Egs. lpO)) 
and (|61l) . the above solution for the vortex core bound state decays as ^ exp(— ^^x^ r) in the long distance regime. On 
the other hand, the contribution from quasiparticles vi^ith the large Fermi momentum to the vortex core state decays 
like ~ exp(— ^r) (note that when /i = 0, the Fermi velocity vp ^ 2A). Thus, for < ft- — A < A, the long distance 
behavior of the vortex core state is dominated by quasiparticles with fc ~ 0, and hence the above approximated 
solution mainly constructed from quasiparticles with A; ~ is valid under this condition. 

In the above derivation of the zero energy Majorana bound state, we have used several approximations. In particular, 
at the stage of matching the solutions for r > Tc and r < rc, we have neglected corrections of order 0{h/{mX^)). 
Actually, such approximations are not essential for the realization of the Majorana zero energy mode, but, rather, 
required by our approximation method of the construction of the zero mode. In fact, the exact solution for the 
zero energy mode in a vortex core should be a superposition of quasiparticles with fc ~ and those with kp ^ Q, 
because the momentum is not a good quantum number in the presence of a vortex core. Since it is quite difficult to 
obtain the exact solution of the zero energy mode constructed from both quasiparticles with fc ^ and those with 
/cf 7^ 0, we approximate it by the bound state mainly formed by quasiparticles with fc ~ 0, neglecting contributions 
of quasiparticles from the large Fermi surface. Because of this approximation, we need the additional approximations 
mentioned above, when we match the solutions for r > Vc and r < Vc- We expect that for the exact solution of 
the zero energy Majorana mode, the weight of quasiparticles from the large Fermi surface may become substantially 
large in the short-distance region in the vicinity of the center of the vortex core, compared to the contributions from 
quasiparticles with ~ 0. If we properly include the mixing with quasiparticles from the large Fermi surface, we may 
be able to match the solutions without such additional approximations. Also, our zero energy solution is not regular 
at r = 0, though it is still normalizable, and physically allowed. We expect that this singular behavior for r ~ is 
also raised by our approximation neglecting the mixing with quasiparticles from the large Fermi surface, which should 
be important for small r. If one takes into account contributions of quasiparticles from the large Fermi surface, it 
may be possible to cure this singular behavior of our solution at r = 0. Nevertheless, when the conditions h <^ mX^ 
and 0</i — A<Ais satisfied, the zero energy solution may be dominated by quasiparticles with fc ~ in the 
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long-distance asymptotic regime sufficiently far away from the center of the vortex core, and our analytical solution 
obtained above may become a good approximation. 

We would like to note that although the zero energy solution given by (|60p and ((6T|) with the asymptotic forms of 
(|D29p and (jD30[) in the Appendix D looks like localized even in the limit of A — > 0, this never means that there is 
a zero energy Majorana bound state in the normal state. In fact, the above zero energy solution is not applicable to 
the normal state without a vortex, because of the following reason. The zero energy solutions in the normal state are 
indeed given by (j58p and (j59p . However, the Hankel function is not regular at r = 0. This singularity is unphysical 
because in the absence of a vortex in the normal state, translation invariance is recovered, and thus, there should not 
be a special point in the coordinate space at which the wave function is singular. This implies that the zero energy 
solution given by ([55)) and ([5^ is not allowed in the normal state without a vortex. Thus, there is no zero energy 
bound state in the limit of A — > 0. The Majorana bound state in a vortex core exists only for hbH^ > A > 0. 

The above analysis may be also extended to the case of d + id- wave pairing straightforwardly, since the d + id-wave 
state, which is the superposition of the drc2_y2-wave state and the dxy-wave state, is an eigen state of the orbital 
angular momentum operator, as in the case of the s-wave pairing, and the treatment for the gap function given above 
is applicable. We obtain the Majorana fermion mode for a nonzero magnetic field /i > in the case of the d -fid-wave 
pairing, which is consistent with the existence of the chiral Majorana edge discussed in Sec lVIl 

We stress again that the Majorana zero-energy mode in a vortex core is formed mainly by the superposition of an 
electron and a hole with the vanishing Fermi momentum fc^r ~ in the long-distance asymptotic regime, provided 
that the energy scale of the SO interaction is sufhciently larger than the Zeeman energy and that 0<ft. — A<A. As 
will be discussed in the following, this property is very important for the stability of the Majorana fermion against 
various sources of decoherence which exist in real materials and may destroy the Majorana fermion acting as a qubit. 

B. Strong stability of the Majorana fermion mode against thermal noise 

For the detection of the non- Abelian anyons and also for the implementation of the topological quantum computation 
utilizing them, it is desirable that the Majorana zero energy state in a vortex core is well separated from excited states, 
the interaction with which may cause decoherence. As was pointed out in ref.^-, in the non- Abelian phase of s-wave 
Rashba superconductors, when the energy scale of the SO interaction is much larger than the Zeeman energy and 
the condition < /iB-ffz — A < A is satisfied, the excitation energy of the vortex core state is of order ^bHz — A, 
which is much larger than the typical size of the excitation energy in the vortex core bound state of weak-coupling 
superconductors, i.e. ^ A'^/Ep^ Thus, the Majorana fermion mode found here is quite stable against thermal noise 
even at moderately low temperatures, to which it is not difficult to access within standard experimental techniques. 
We, here, explain the origin of the strong stability of the Majorana mode in more details. The excitation energy in the 
vortex core is due to the kinetic energy of quasiparticles in the Andreev bound state, which stems from the derivative 
term in Eas. (jDl[ ) and (|D2I) in the Appendix D. We restrict the following argument within the case that the magnitude 
of the SO interaction is much larger than the Zeeman energy, and that the condition < hbH^ — A < A is satisfied. 
In this case, the Majorana solution is constructed mainly from quasiparticles with fc ^ 0, as clarified in the previous 
section. Then, the first order derivative terms 2A(^ ± give leading contributions to the kinetic energy. On the 
other hand, from the solution of the BdG equation (feot and (|6T|) . we see that the characteristic size of the vortex core 
is ^corc ^ 2A/(/iB^^2 ^ A). Thus, the excitation energy is of the order ^ 2A/ (2A/ {hbHz — A)) ^ /is^^z — A. This large 
magnitude of the excitation energy implies that the Majorana zero energy mode in the s-wave Rashba superconductor 
with ^.bHz > a is significantly stable against thermal noise, compared to chiral p + ip superconductors. Also, such 
large excitation energy ensures that the experimental detection of the non- Abelian anyons is quite feasible for our 
system. The origin of the strong stability of the Majorana fermion mode is deeply related to the fact that it is mainly 
constructed from quasiparticles in the Dirac cone at the F point in the Brillouin zone for < hbHz — A < A in the 
long-distance asymptotic regime, as mentioned above. Since the Dirac cone has a vanishing Fermi momentum, the 
kinetic energy is dominated by the SO interaction which is of order A/^corc rather than the standard kinetic energy 
term of order {2m). This feature leads to the strong stability of the Majorana zero energy mode. It is noted that 
the robustness of the Majorana fermion mode in the s-wave Rashba superconductor was also pointed out in refi^^ 
from a different point of view. 

C. Stability against decoherence due to inter- vortex tunneling 

It has been proposed that the Majorana fermion modes in superconductors can be utilized as decoherence-free qubits, 
which enable us the construction of the fault-tolerant topological quantum computer Two Majorana fermions, 

say 71 and 72, constitute one complex fermion state described by '0 = 71 + 172, which is occupied or unoccupied. This 
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doubly-degenerate state stores a qubit non-locally, which is protected against any local perturbations, as long as the 
distance between two vortices, each of which contains one Majorana mode is sufficiently large. One crucial obstruction 
to this scheme is the decoherence raised by inter-vortex tunneling; tunneling processes between two Majorana modes 
in two vortices lift the degeneracy, leading to decoherence. In particular, it was pointed out by Cheng et ali^ that 
the energy of the complex fermion -0 exhibits quantum oscillation as a function of the spatial separation r between 
two vortices as ^ cos kpr with kp the Fermi momentum, and thus takes both positive and negative values depending 
on r. This rapid change of the sign of the energy seriously flaws the initialization and the readout of the qubit. 

Here, we discuss the drastic suppression of the quantum oscillation for inter- vortex tunneling in the Rashba s-wave 
superconductors in a particular parameter region. As explained in Sec lVII M the zero energy Majorana bound state 
in a vortex core obtained above consists of two contributions; one from quasiparticles with A: ~ and the other from 
quasiparticles with the large kp. The former contribution has the characteristic length scale of order A/(/j,b^^z ^ 
as clarified by the solution (pO]) and (pT|) . On the other hand, for the latter contribution, the characteristic length scale 
is of order A/ A, since the Fermi velocity is of order A when ^ ~ in Eg. ((57)) . Then, in the case of < /ib-H^z ~ A < A, 
the Majorana zero mode is mainly formed by quasiparticles with the vanishing Fermi momentum hp ^ in the 
long-distance asymptotic regime. That is, in the tunneling process between two vortices separated by the distance R, 
the overlap of the component with A: ~ 0, which is of order exp(— i?(/iB^z — ^)/-^) is much larger in the magnitude 
than the oscillating contribution from the large Fermi surface which is of order exp(— i?A/A) cos(fciri?) for large R. 
The inter- vortex tunneling mediated via quasiparticles with kp ^ does not involve the quantum oscillation, making 
a sharp contrast to Majorana modes found in p + ip superconductors. As a result, the energy of a complex fermion 
made of the two Majorana fermions is dominated by the non-oscillating part from quasiparticles with /c ~ 0, and thus, 
the decoherence due to the quantum oscillation is suppressed. It is noted that this protection mechanism also works 
for a Majorana fermion mode in the proximity between a topological insulator and an s-wave superconductor,^'^ as 
long as the chemical potential is properly tuned to realize kp ^ for the surface Dirac cone. Since momentum is not 
conserved in the vicinity of a vortex, there may be hybridization between the Majorana state in a vortex core and 
quasiparticles with the finite Fermi momentum, which raises inter- vortex tunneling involving the quantum oscillation. 
However, the analysis of the zero energy Majorana mode in a vortex core in the previous subsections implies that we 
can construct the zero energy Majorana state which is mainly formed by quasiparticles with k ^ 0, provided that the 
energy scale of the SO interaction is much larger than the Zeeman energy, and that < (J-bHz — A < A. Thus, it may 
be possible to realize the Majorana fermion in the vortex core of the Rashba superconductor which is stable against 
decoherence due to the inter- vortex tunneling. 

VIII. TOPOLOGICAL DENSITY WAVE STATES AND CHARGE FRACTIONALIZATION 

The above argument for the non-Abelian topological order in s-wave superconductors implies that topological order 
is realizable in conventional spin (or charge) density wave states. The Hamiltonian for the s-wave superconductivity 
on a bipartite lattice is mapped to the Hamiltonian for the spin density wave (SDW) state or the charge density wave 
(CDW) state by changing the basis of fermion fields. Thus, an s-wave superconducting state with the topological 
order can be mapped to a density wave state with a certain topological order. As shown below, the topological order in 
the density wave state is Abelian, and quasiparticles in this topological order are not Majorana fermions, but possess 
U(l) charge. However, the quasiparticles exhibit charge fractionalization which characterizes the Abelian topological 
order. 



We, first, consider the SDW state with the order parameter As = {cl^Ck+Q^^) — — (c^c/j+qj^), where Q is the 
ordering wave number vector. We also assume that there are the Rashba-type SO interaction and the Zeeman 
magnetic field h = ^i^Hz- Then, the mean field Hamiltonian is given by 



A. Topological spin density wave state 
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TABLE IV: The TKNN integer in the magnetic BZ, JmTKNN, and the winding number I{ky) for 2D SDW state 



and ^'1" = (cfe^, Cfe^, Cfc+Q^, Cfc+Q^). In the following, we assume the perfect nesting condition of the energy band 
e{k + Q) = —e{k). This situation is realized in the case of the half-filling electron density for our model on the 
square lattice with the nearest- neighbor hopping. The nesting vector is Q = (±7r,±7r). Furthermore, we postulate 
^o{k + Q) = Co{k). This condition is satisfied when Co{k) — (sin2fcj,, — sin2fca;) in the case Q — (±7r, ±7r). Then, 
the Hamiltonian (j63p has the same form as that of the Rashba s-wave superconductor considered in the previous 
sections, and topological order is realized when the condition h > Ag is fulfilled. 

In a similar manner as Sec|Vl topological numbers of the SDW state are calculated. The energy gap of the system 
closes one of the following conditions are satisfied, 

AI = {pbH,)\ At^ + Al = {pBH,f, 16^2 + Ai = (A^Bi?.)', (64) 
and the winding number I{ky) is calculated as 

I{ky) = \ [sgn[e(-7r/2, k^f - + A^] - sgn[£(0, k^f ~ + A^] 

+sgn[£(^/2, kyf - (MBi?.)' + A2] - sgn[£(7r, kyf - (fiBH.f + A^]] . (65) 

Because T~1-sdw(~^) — ^SDw(fc) at ky — 0, ±7r/2,7r, the winding number I{ky) is defined at these values of ky. We 
summarize the winding number of the SDW state in Table IIVI The TKNN integer in the magnetic Brillouin zone, 
which we denote as /mTKNN, is also given in Table ITVl 

In a similar manner as s-wave NCS superconductors, gapless edge modes appear in the NCS SDW state considered 
above when h > As- See FiglS] However, in comparison with the case of s-wave superconductivity, we should taken 
into account only half of the zero energy modes, since the Brillouin zone is folded to the magnetic Brillouin zone in the 
SDW state, and as a result, the fc-point k = (tt, 0) is equivalent to fc = (0, tt). In other words, when (— l)^>nTKNN _ 
we have odd numbers of zero energy edge modes in the SDW state. In FiglHl we see three zero energy modes. The 
existence of odd numbers of zero energy edge modes implies that in the vortex core of the SDW order Age'^, there 
are odd numbers of zero energy modes, because of the bulk-edge correspondence. These zero energy modes have U(l) 
charge. As in the case of the Su-Schrieffer-Heager model^^ these isolated zero energy modes with U(l) charge lead 
to the charge fractionalization4i That is, when these three zero modes are occupied by electrons, the charge carried 
by the vortex core is Q = 3e/2. 



B. Topological charge density wave state 

The above consideration for the topological SDW state is also applicable to the CDW state. We consider the CDW 
state with the order parameter Ac = (cj^^Cfc+g^) = {cj^^Ck+Qi) in the case with the Rashba SO interaction. The mean 
field Hamiltonian for the Rashba CDW state is 

HcDW ^IY. *k^CDw(fc)«'cfc, (66) 

k 



with 



HcDw(fc) - ( '^''^ ' + "^''^^^ ■ " '^^^^ 1 (67) 

I -iAc<Jy e{k + Q) + ha, + aCo{k + Q) ■ <J* I 
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FIG. 9: The energy spectra of the NCS SDW state with edges at i^: = and = 30. Here ky denotes the momentum in 
the {/-direction, and ky G [— 7r,7r]. We take i;=l = a = l = As = 1. The Zeeman magnetic field Hz is fiaHz = f.5. At 
ky = 0, ±7r/2, TT, we have gapless states on each edge. These gapless edge modes are chiral. 



and ^'^j, — (cfc^, Cfej., — Cfe+Qj., Cfe+Q-t-). We assume the prefect nesting condition e{k + Q) = —£{k) again, and also 
assume that Coi^k + Q) = — £o(fc)- The condition for Co{k) is satisfied when Co{k) — (sin ky, — sin k^) in the case of 
Q = (±7r, ±7r). Eq. (|67l) is formally equivalent to the Hamiltonian of the Rashba s-wave superconductor. In this CDW 
state, the topological ordered phase with a single gapless edge mode realizes. We depict an example of the energy 
spectra for the system with open boundary edges in the case of ft, > Ac in Fig llOl There are two zero energy edge 
modes at fcy = and tt. Since the Brillouin zone is folded by the CDW order with the ordering vector Q = (tt, tt), the 
zero mode at ky = tt is equivalent to that at fcj, = 0. Thus there is only one zero mode. In this Abelian topological 
phase, the charge fractionalization can occur, as in the case of the NCS SDW state. When there is a vortex of the 
CDW order, the zero mode of which is occupied by one electron, the vortex carries the fractional charge e/2. 

IX. DISCUSSION AND CONCLUSION 
A. Realization scheme for the non-Abelian topological order in spin-singlet superconductors 

In this section, we discuss possible realization schemes for the non-Abelian topological order considered in the 
previous sections. Recently, it was discussed by several authors that s-wave Rashba superconductors with the condition 
HbHz > A are realizable in heterostructure semiconductor devices i^^'^^ We, here, present two more possible schemes; 
one based on quasi-two-dimensional bulk superconductors, and the other utilizing ultracold atoms. 

1. Possible realization in strongly correlated electron systems 

The most crucial hurdle for the realization of the non-Abelian topological order in the Rashba superconductors 
is to stabilize the superconducting state under applied strong magnetic fields ^bHz > A. As mentioned before, 
the Pauli depairing effect due to magnetic fields is absent when the SO splitting of the Fermi surface is sufficiently 
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FIG. 10: The energy spectra of the NCS CDW state with edges at i^, = and = 30. Here ky denotes the momentum in the 
y-direction, and ky £ [— 7r,7r]. We take t = 1, = 0, A = 0.5, and Ac = 1. The Zeeman magnetic field H^, is /xb^^z = 1.5. 



larger than the superconducting gap and the Zeeman energy .^'^'^^ Such large SO interactions are typically realized 
in noncentrosymmetric materials. Then, the issue which should be addressed here is how to suppress the orbital 
depairing effect of magnetic fields. The typical size of the orbital limiting field is given by ^bH°2^ ~ /S.^ /{zEp). 
Here z is the mass normalization factor. If the mass enhancement factor 1/z is sufficiently large, it is possible to 
attain ji-gH"^ > A. For instance, for strongly correlated electron systems such heavy fcrmion superconductors, 1/z 
can reach to be ^ 100. Thus, it is not difficult to fulfill the condition ^-qH"^ > A even when A/Ep ^ 0.01. The 
above argument implies that noncentrosymmetric heavy fermion superconductors may be good candidates for the 
realization of the non-Abelian topological order. Unfortunately, to this date, there are no noncentrosymmetric heavy 
fermion superconductors with no gap-nodes. Experimental studies suggest that all heavy fermion superconductors 
without inversion symmetry discovered so far such as CePtsSi and CeRh(Iir)Si3 possess nodes of the superconducting 
gap^^-il excitations from which may destabilize the topological order. 

2. Possible realization in ultracold fermionic atoms 

We now discuss an experimental scheme for the realization of the non-Abelian topological order in ultracold fermionic 
atoms. We can use the Feshbach resonance in the s-wave channel for the formation of the s-wave Cooper pairs in this 
system^i^. For superfluid states of charge neutral atoms, there is no orbital depairing effect due to applied magnetic 
fields. Moreover, it was recently proposed by several authors that effective SO interactions acting on charge neutral 
atoms can be generated by spatially varying laser fields^"—. When the effective SO interaction is the Rashba type 
(or more generally, anti-symmetric), and its magnitude is sufficiently larger than the Zeeman energy due to an applied 
magnetic field, the Pauli depairing effect is suppressed, as mentioned before. Thus, the condition for the non-Abelian 
topological order fisHz > A can be easily fulfilled. As discussed in the previous subsections, the non-Abelian anyons 
are stable for sufficiently low energies <^ min{/XB-ffz — A, A}. Since the BCS gap A can be tuned to be large, i.e. 
A ~ Ep, by using the s-wave Feshbach resonance, the realization of the non-Abelian anyons in this scheme is quite 
promising. In the following, we present two schemes for the realization of the topological phase in ultracold fermionic 
atoms with laser generated SO interaction. 

The first scheme utilizes an optical lattice with laser-assisted tunneling of fermionic atoms. We consider fermionic 
atoms loaded in a two-dimensional (2D) periodic optical lattice. The atoms occupy doubly-degenerate Zeeman 
levels of the hyperfine ground state manifolds, which correspond to the spin up and spin down states of electrons. We 
introduce the Zeeman field to lift the degeneracy. We assume that standard tunneling of atoms between sites due to 
kinetic energy is suppressed by the large depth of the optical lattice potential. Tunneling of atoms between neighboring 
sites along the i/-direction (v = x, y) which conserves spins is caused by laser beams via optical Raman transitionsi^i^ 
In addition, tunneling which accompanies spin flip is also driven by two Raman lasers^SS The Rabi frequency of the 
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laser fi^i (^^2) is resonant for transition from the spin- up state to the spin-down state for the tunnehng between 
neighboring sites in the forward (backward) z/-direction. Furthermore, the confining optical potential is tilted along 
both the a;-direction and the y-direction to assure that the forward and backward tunneling processes are respectively 
induced by the lasers with the different Rabi frequencies n^i and ^1^2- The tilting potential for the a:-direction A^; 
is different from that for y-direction Ay to ensure that tunneling accompanying spin flip along the ?/(a;)-direction is 
not raised by the lasers with ^x{y)i.2- It is also assumed that the detuning from excited states for optical Raman 
transitions is much larger than A^(j^), and thus the spatial variation of the amplitudes of the Rabi frequencies due 
to the tilting potential is negligible. To realize the Rashba SO interaction for the two Zeeman levels, we choose the 
phases of the lasers as follows. The lasers are propagating along the z-direction with an oscillating factor e*'^''^. The 
Rabi frequency ^1^2 is expressed as ^13,2 — \rtx2\s^^'^- The phase of the laser ft^i is shifted by n from that of ^1x2, 
and r2a;2 — —^xi holds. The phase of flyi (f2j,2) is shifted by — 7r/2 (7r/2), and ily2 — —iflxi, ^y2 — ^^yi- Then, the 
laser-induced tunneling term which accompanies spin flip is expressed by 



with v = X, y, and Cx — I and Cy = —i. Since we consider the 2D xy-plane with z = 0, ai, is real. For ax — cty, 
Eg. (155)) is the Rashba SO interaction. 

In the second scheme, we employ the idea that an effective SO interaction is created by utilizing the dark states 
generated by spatially-varying laser fields, as proposed in refsi^i^ We can introduce a vortex of the superfiuid order 
parameter by changing the topology of the shape of a trapping potential. Since the underlying topological effective 
field theory is the SU(2)2 Chcrn-Simons theory ; '^'^'^^ a vortex with a zero energy Majorana mode is equivalent to a hole 
pierced in the system provided that the total number of the holes in the system are odd. The trapping potential with 
holes can be prepared by using a holographically engineered laser technique. Furthermore, the dynamical motions 
of holes are possible through a computer-generated hologram; i.e. the positions of holes can be moved by changing 
temporally the shape of the confining potential which may be achieved by preprogrammed hologram. Remarkably, 
this enables the braiding of vortices (holes), which amounts to a quantum gate for fault-tolerant quantum computation 
based upon the manipulation of non-Abelian anyons-i2.'«i2. 



We have verified that the non-Abelian topological order, which is characterized by the existence of chiral Majorana 
edge modes and a Majorana fermion in a vortex core, can be realized in almost all classes of fully-gapped spin- 
singlet superconductors with anti-symmetric SO interactions such as the Rashba SO interaction in the case with 
the Zeeman magnetic field. Majorana fermions in superconductors have been, recently, attracting much attention 
in connection with the topological quantum computation. The idea of the topological quantum computation has 
been examined for the non-Abelian fractional quantum Hall state, "'^^'^'^ spin-triplet chiral p-wave superconductors,— 
and the proximity between a topological insulator and a conventional superconductor. ^^'■^^ Our results in the present 
paper provide another category of promising systems for the realization of the topological quantum computation. 
We would like to stress that our systems in the normal state are topologically trivial; i.e. the combination of the 
conventional superconducting order and the conventional anti-symmetric SO interaction with the Zeeman field gives 
rise to a highly nontrivial topological phase. In our scenario, the existence of the Zeeman energy, the magnitude of 
which exceeds the size of the superconducting gap opened at time-reversal invariant /c-points in the Brillouin zone 
where the SO interaction vanishes, is important. At these fc-points, there are Dirac cones in the absence of magnetic 
fields. When the magnetic field is switched on, the effective gap of electrons in the vicinity of the Dirac cone is 
given by A — ^bHz- As the magnetic field increases, and ^bHz exceeds A, the sign of the effective gap changes, 
and topological phase transition occurs. As indicated by the analysis of Majorana fermion modes in open boundary 
edges and in a vortex core in Sec lVIl and Sec lVIIl and in the topological argument based on the winding number in 
Sec|Vl quasiparticles in the vicinity of the Dirac cone, i.e. time-reversal invariant /c-points, play an important role 
for the realization of the topological quantum phase transition. In the case of the s-wave pairing state, the large 
magnetic field satisfying (IbHz > A is required to realize the non-Abelian topological order. For conventional bulk 
weak-coupling superconductors, it is difficult to fulfill this condition, though there are some proposals to realize this 
system in heterostructure devices, for which ^jlbHz > A is attainable i^^'^^ We have proposed promising schemes for 
the realization of the condition (IbHz > A in the bulk s-wave pairing state; one based on strongly correlated electron 





(68) 



B. Conclusion 
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systems, and the other utihzing ultracold fermionic atoms. In contrast to the s-wave pairing state, in the case of the 
d + id wave pairing state, a smah magnetic field larger than the lower critical field is sufficient to fulfill the condition 
IJbHz > Afe in the vicinity of time-reversal invariant fc-points, and thus the realization of the non-Abelian topological 
phase is easier than the s-wave case. However, unfortunately, up to our knowledge, there is no superconductor for 
which the d + ic?-wave pairing state is experimentally established, though it was suggested by some authors that this 
pairing state may be realized in sodium cobalt oxide superconductors Naa;Co02-yH20 i^^i^^ It is an interesting open 
issue to pursue the experimental detection of Majorana fermions in such superconductors mentioned above. We have 
also clarified that, in certain parameter regions, the Majorana fermion realized in Rashba spin-singlet superconductors 
is stable against various sources of decoherence such as thermal fluctuations which exist in real systems generically. 
Because of this feature, the Rashba spin-singlet superconductor may be a promising candidate for the realization of 
the topological quantum computation. 

As a byproduct of our analysis, we have also found that an Abelian topological order, which supports the existence 
of gapless Dirac fermions in edges and the charge fractionalization, is realizable in the conventional SDW or CDW 
state with an anti-symmetric SO interaction. 
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Appendix A: TKNN number and the dual transformation 

In this appendix, we show that the dual transformation (fTTj) does not change the TKNN number, and hence the 
topological properties of the original Hamiltonian (jlj and the dual Hamiltonian are the same. Let us consider 
the BdG equations for the original Hamiltonian H{k) and that for the dual one T-Luik), 

n{k)\Mk)} = E4k)\Mk)), nu{kM^{k)) = E,,{k)\<f>^{k)). (ai) 

Since the dual transformation is a unitary transformation, the energy spectrum of the dual Hamiltonian is the same 
as that of the original one. Moreover, the eigen state |0j^(A:)) for the dual Hamiltonian is related to the eigen state 
\(pn{k)) for the original one as 

\cjy^ik))=D\Mk)), (A2) 

with D in ()I2|) . As _D is a constant unitary matrix, the gauge field constructed from the occupied states for the dual 
Hamiltonian is the same as that for the original one, 

4"^°(fc)=* E {<Pn(.k)\^,,^^ik))^^ E {Mk)\dMk)) ^ A'r\k). (A3) 

-E„(fc)<0 E„{k)<0 

This relation immediately means that the TKNN number (jl5p is invariant under the dual transformation. 



Appendix B: TKNN number and winding number 

1. TKNN number for TRB topological superconductors 

In this appendix, we summarize useful properties of the TKNN number for TRB topological superconductors.— 
Let us first consider the BdG equation, 

H(fc)|0„(fc)) = £;„(fc)|0„(fc)), (Bi) 

where the BdG Hamiltonian 'H(fc) has the particle-hole symmetry, 

rn{k)Ti = -n{-k)*, (B2) 
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with the 4x4 matrix T 



r = 



l2x2 
l2x2 



(B3) 



From the particle-hole symmetry, we can say that if |0„(/c)) is a positive energy state with En{k) > 0, then r|(/>* (— fc)) 
is a negative energy state with —En{—k). Therefore, in the following, we use a positive (negative) n to represent a 
positive (negative) energy state, and set 



|(/._„(fc)) = r|0:(-fc)). 

For the ground state in a superconductor, the negative energy states are occupied. 
Now we introduce the following "gauge fields" A'f^^k), 



From (|B4p . the gauge fields A^j^\k) satisfy 



4+\k) ^ 4-\-k). 



(B4) 



(B5) 



(B6) 



It is also found that their sum Ai{k) = A'l^\k) + A'i \k) is given by a total derivative of a function. To see this, we 
rewrite 0„(fc) in the components. 



and introduce the 4x4 unitary matrix W{k) as 



/ 4'Uk) \ 
ct^lik) 

\ ^iik) J 



(B7) 



WaAk)=rn{k). 

Then Aiik) is rewritten as the total derivative of iln[detiy(fc)], 

A,{k) = it! [w\k)dk,W{k)\ = iafedn[detW^(fc)]. 
This equation implies that the field strength of Ai{k) is identically zero, 

F{k)=e'^dk,Aj{k)^{). 

Combining this with (jB6[) . we find that the field strength F'^^\k) of A^^\k) satisfies 

J-(±)(fc) =^(±)(-fc). 
By using the field strength of the occupied state, the TKNN number is define as 

1 



'TKNN 



2ti 



rp2 



dPkF^-\k), 



(B8) 
(B9) 
(BIO) 

(Bll) 

(B12) 



where is the first Brillouin zone in the momentum space. Using the relation (jBlip and the Stokes' theorem, we 
obtain 



'TKNN 



<fkF^-\k) 
dk,Ai-\k,,Q) - I dfc,AW(fc„7r) 



(B13) 



where is the upper half plane of T^. As is shown in Appendix IB 3[ this formula enables us to connect the Chern 
number to the winding number defined in (|2ip . 
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2. winding number 

As was discussed in Sec|Vl the winding number I{ky) is defined by taking the basis where the BdG Hamihonian 
has the following particular form^ 

nik)=( , '^^''A. (B14) 



gt(fc) 



By using q{k), the winding number is defined as 



I{ky) = ^ r rf^.tr [q-Hk)dkMk) - q^'\k)dkj (k)] , (B15) 



which is equivalently rewritten as 



I{ky) - / dk,tr[q{k)dk,q-\k)] 



/ dA:,afe.ln[det9(fc)]. (B16) 



Now we derive a useful formula to evaluate the winding number I{ky). Denote the real and imaginary parts of 
det5(fc) as TJiiik) and m2(fc), respectively, 



Then, I{k) is rewritten as 



where rhi{k) is given by 



dctq(k) = mi{k) + im2{k). (B17) 



Hk) = ^ I dk,t''m,{k)dk^rh,{k), (B18) 



y/mi[ky + 777,2 (fcj 



To evaluate (jB18l) . we use the technique developed in^S,. From the topological nature of I{ky), we can rescale one 
of mi{k)s, say 77ii(fc), as mi{k) — ?• ami{k) (a <C 1) without changing the value of I{ky). Then it is found that only 
neighborhoods of k^ satisfying 7772 (fc^, ky) — contribute to I{ky) if a is small enough. By expanding rai{k) as 

777l(fc) =7r7l(fc°,fcj,) + --- , 7772 (fc) =9fe^7772(fcS,fcj,)(fc:r-A:°) + --- , (B20) 

the contribution from /c^ is calculated as 

isgn[777l(fc°,fcj,)] •Sgn[afc^7772(fc°,fcj,)]. (B21) 

Summing up the contribution of all zeros, we obtain 

I{ky)^ E \sgn[m^{kl,ky)]-sgn[dk^m2{kl,ky)]. (B22) 

Exchanging mi{k) with 7772 (fc) in the above argument, we also have 

Iiky) = - \^Sn[9k^mi{k°,ky)].-sgn[m2{kl,ky)]. (B23) 

{kO:mi(kO,ky)=0} 
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3. The relation between the TKNN number and the winding number 

Here we prove the relation (|22l) between the TKNN number and the winding number. As was shown in (|B13[) . the 
TKNN number for a TRB superconductor is evaluated by the line integral 



2^ rTT ^ ^ 



(B24) 

with ky = 01 ky — TT. From the particle-hole symmetry, this line integral itself is a Z2 topological numberi ^^i^° In 
this section, we relate this line integral to the winding number defined in the previous subsection. 
First, consider the eigen equation 

q{k)q^{k)\ur,{k)) = En{kf\u„{k)), (B25) 
where q{k) is given by an off-diagonal component of the Hamiltonian (|B14[) . Using the components of |u„(/c)), 



(B26) 



we define the 2x2 unitary matrix U{k), 

Uik)a,n^Kik). (B27) 

Then the eigen equation is recast into 

q{k)q^k)U{k) ^U{k)A{k), (B28) 

where A(A:) is given by 

A(fc) - diag{Ei{kf,E2ik)^). (B29) 

Equation (|B28| yields 

g-i(fc) = q^k)U{k)A-\k)UHk). (B30) 
By using q{k) and |u„(fc)) in the above, the occupied state \(j)^^\k)) of the Hamiltonian (|B14[) is given by 

u(-)(fc)) = ^ ( , '"""(^^^ I , (B31) 

^ V2\qHk)\u,,{k))/Er.{k) y ' ^ ' 

with negative n. Here \un{k)) is normalized as (u„(/c)|u„(A:)) = 1. Thus the gauge field \k) is calculated as 

4~\k) = ^J2{<pi-Hk)\^,A-\k)) 

n<0 

= i 

n<0 n<0 ri<0 ' 

n<0 ' 



Y^{u^{k)\du^{k)}+iY,^^\jT^<*{k)qab{k)dk^^^^^^ (-i- 

n<0 n<0 ' ri<0 ^ "W, 

itr [u\k)dk,U{k))] +i^tT [q{k)dk,q^U{k)A-\k)U\k)] + iY,lE^ik)dk, f-^Tn) 

itr [[/t(fc)9fc^[/(fc))] + iitr [g(fc)afe.g-i(fc)] + 4 E ^fe.lni?„(fc), (B32) 



2 

n<0 



where we have used (jB30l) in the last line of the above equation. We also notice here that the unitary transformation 
used to obtain the Hamiltonian in the form (IB14p is accomplished by a constant unitary matrix, so it does not change 

the value of A\ \k). Substituting this into (IB13|) . we obtain 



/tknn = /(O) - /(tt) -I- - / dkxir[u\k)dk^U{k)\ -- / dk.Ar [u\k)dk^U [k)] 



(B33) 



Since we have 



dkxiy[u\k)dk^U{k)\=^ j dk^dklndeiU {k) ^ 2N , (B34) 
with an integer iV, the last two terms in the right-hand side of (jB33l) are even integers. Therefore, we obtain 
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Appendix C: TKNN number and chirality basis 

In Sec|Vl it is discussed that the origin of the non-Abehan topological order in spin-singlet Rashba superconductors 
is understood in terms of mapping to spinless odd-parity superconductors which is derived from the chirality basis 
representation. This mapping was first considered in ref.— , though the parameter region in which this mapping is 
applicable was not fully elucidated in refi^. Actually, this mapping can be used when the Zeeman energy is much 
larger than the superconducting gap. In this appendix, we will show that this mapping does not change the topological 
number of the original system. This property is not trivial, since the mapping depends on wave numbers and hence 
the band structure of electrons in a non-trivial way. 



1. chirality basis 

For simplicity, we suppose that /ib^z > in the following. In our Hamiltonian Q, the normal dispersion of electron 
is determined by 

£{k)^e{k)~fiBH,a,+aCo{k)-<J. (CI) 
In the chirality basis, it is diagonalized as 

m = U{k)(^'+ll'^ ^ \f^^^uHk), (C2) 

where e±{k) — e{k) ± Ae(fc), and Ae(fc) is given by 

^e{k) = ^{aCoikW + {pbH,)\ (C3) 

and U{k) is 

jji^A^ 1 ( aCoxik) - iaCoy{k) -Ae(fc) - ^b-^z \ 

^2Ae{k){Ae{k) + fiBH,) \ Ae(fc) + ^Bi?. aCo^ik) + iaCoy{k) J ' 

From the following unitary transformation, 

n{k)^G{kyn{k)G{k), (C5) 

with 



/ t/t(fc) 
1^ U^{-k) 



Gik)^ ' ,1, (C6) 



it is found that the BdG Hamiltonian T-L{k) in the chirality basis is given by 

^ / e(fc) + A£(fc)a, zA(fc) [[/t(fc)a,t/*(-fc)] \ 

^' \-iA*{k)[uT{-k)ayU{k)\ -e{k) - Ae{k)a, J' ' 

Therefore, the gap function Aa-a' (k) in the chirality basis becomes 



A^.'ik) 

tAik) [uHk)ayU*{-k)]^^, 



(aCoxik) + iaCoy{k))A{k) ^Bi?zA(fc) 



Ae(fc) \ -^Bi?.A(fc) (aCoAk) - taCoy{k))A{k) 



(C8) 



This equation indicates that when the original gap function A(fc) is even-parity, the odd-parity gap functions (the 
diagonal terms of (jCSP ) are induced in the chirality basis due to the existence of the SO interaction. However, this 
never means that the topological class is always the same as that of spin-triplet superconductors, because the off- 
diagonal terms of (jCSp corresponding to the inter-band pairing may change its topological property. As a matter of 
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fact, the non-Abelian topological order appears when one of the two bands s±{k) is gapped by the large Zeeman field 
satisfying (J-bH^ S> |A(A;)|, and only one band survives in the low energy region. More precisely, in the parameter 
regions of the non-Abelian phase shown in Tables I, II, and III in Sec|Vl when /i < 0, the Zeeman term generates a 
large gap in the band e+(fc), leaving only one band £-(fc) in the low energy region, and hence a spinless odd-parity 
superconductor is realized for the band e_(fc). The effective Hamiltonian for the spinless chiral odd-parity 

superconductor is obtained by integrating out fermion fields for the high energy massive band: 

n(k)=( {aCoAk)-^aCoy{mA{k)/Ae{k))\ 
\{aCoAk) + iaJ^oy{m^*{k)/Ae{k)) -e_(fc) J' 

In a similar manner, in the non-Abelian phase with fi > 0, the Zeeman term generates a large gap in the band £_(fc), 
leaving only one band e+(fc) in the low energy region, and thus, the following spinless odd-parity superconductor is 
realized for the band £+(fc), 

~ / e+{k) ia£oAk)+iaCoy{k)){A{k)/Ae{k))\ 

+ ' ' \{aC,,{k)~ta£oy{k)){A*{k)/Ae{k)) ~e+{k) )' ^ ' 

In both cases, the Rashba s-wave (d -\- id-wave) superconductor in this situation is mapped into spinless chiral p-wave 
(/-wave or p-wave) superconductor. 

The above consideration based on the chirality basis is also useful for understanding the origin of the Abelian 
topological order presented in Sec lVIIll The Abelian order realizes in the vicinity of the half-filling /i « 0. In this 
case, there are one particle-like band and one hole-like band. In addition, there are Dirac cones at A; = (±7r, 0), 
(0, ±7r). When hbUz ^ 1^(^)1, the Dirac cone bands have a large gap ~ HbHz, and can be integrated out in the 
low energy region. Then, there remain one particle-like band and one hole-like band. Furthermore, when a|£o(^)| ^ 
HbHz ^ |A(fc)|, the inter-band pairs (off-diagonal terms of Eq. (jC8[) ) are negligibly small compared to the intra-band 
pairs. Since the chiral gapless edge state associated with H+{k) which corresponds to the particle-like band and 
that with H^{k) corresponding to the hole-like band have the same chirality, propagating in the same direction, the 
perturbation due to the inter-band pairs does not raise gap in the two edge modes. This implies that the system is 
mapped to two decoupled spinless chiral odd-parity superconductors for ck|£o(fc)| 3> /is-ffz 3> |A(/c)|. 

By contrast, in the parameter regions where there is no topological order, there are two particle-like bands (/i < 0) or 
two hole- like bands (/i > 0). In such cases, the edge modes associated with i?+(fc) and H-{k) propagate, respectively, 
in the opposite directions. Thus, even a small perturbation due to the inter-band pairing terms of (jC8| gives rise to 
a gap in the edge excitations, and hence there is no topological order. 



2. TKNN number in the chirality basis 

Here we will show that the unitary transformation (|C5|) for the chirality basis does not change the TKNN number 
in the presence of the Zeeman magnetic field, and hence the topological properties of the original Hamiltonian 'H{k) 
and the Hamiltonian H(fc) in the chirality basis are the same. 

Let us first rewrite the formula (jB13l) in a more convenient form. Using (jB6|) . the TKNN number is recast into 

1 r 

/tknn = -/ dk^[A^{k^,{))- A^{k^,Ti)]. (Cll) 

Jo 

We will use this formula to prove the above statement. 

As was shown in Sec IC 1[ H{k) is related to the the original Hamiltonian H{k) as follows, 

n{k) = G{kyn{k)G{k) (C12) 
Therefore, the eigenstate |^„(/c)) for T-L{k) is also related to the eigenstate \(f>n{k)) for T-L{k) as 

\Mk)) ^G{k)\^Jk)). (C13) 
When I-IbHz 0, G{k) is non-singular, thus we have 

Mk) = ^Y.^<l>n(k)\dMk)) 

n 

= iJ2{Mk)\dk.Mk)) - »tr [G^dkMk)] 

n 

= A,{k) - idkM [detG(A:)] = Ai{k), (C14) 
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where we have used detG'(fc) = 1. Therefore, from (IC11|) . it is found that the TKNN number remains the same in the 
chirahty basis. 



Appendix D: An approximated solution for the Majorana zero energy mode in a vortex core 



We, here, present a derivation of an approximated solution for the Majorana zero energy mode in a vortex core 
considered in Sec I VIII Zero energy solutions of the BdG equation generally satisfy the condition = w* , because of 
the particle-hole symmetry of the BdG Hamiltonian. Thus, the BdG equation ([SS]) for i? = can be recast into the 
following two equations for and ui. 



0)2 



2m V dr'^ r dr 86 Ar'^ ^ 



1 0)2 



d 



A^iJ = (Dl) 



1/92 19 1 92 n d n ^ 



- 2Ae^'' I a; + ; - ^ ) - = (D2) 



Here A = for r < r^, and A 7^ for r > r^, with r^, ^ ^. 

We first consider the solution for r < r^. Then, we put A = in Egs. ljDip and (ID2I) . We examine the following 
form of the solution, 



ut(r,0) =e-^^e^('5-i)^/<(r), Mr,9) 



(D3) 



Here /3 is a constant which will be determined later. Substituting (|D3I) into (jDl[) and (jD2|) . we have the equations for 
/< 



t(4.)' 



2m \ 8r'^ r 8r ' ' 



(D4) 



2m I r 8r ' 



(D5) 



We search for the solutions of (|D4[) and (jPSP in the form f^{r) = A^Zp_i[ar) and f^{r) = A^Zp{ar), where Z^{ar) 
is the Bessel function and a is a constant. Substituting these expressions into Eas. (|D4p and (|D5|) and using the 
following relations for the Bessel functions Zy(ar), 



8Z^{ar) ^ , \ r \ ^ r \ 

= -Zy{ar) - aZ^+i{ar) = aZy^i{ar) Z^{ar), 

or r r 



(D6) 



we obtain. 



^ 8P_ _^ ld_ 

Qj,2 J, Qj, 



+ 2m(/Lt + h) — 4mAa 



Zp^i(ar) = 0, 



(D7) 



1 8 p 



7— r + - ^ 7 + 2m(^ -h) - -imXa— 



A, 



Zfi{ar)=Q. 



(D8) 



Eas. (jD7l ) and (jPSp are actually the Bessel differential equations with the solutions Z^_i(Q:r) and Zf^{ar), respectively, 
provided that the following relation is satisfied. 



A I Af 
2m(/i + h) — 4mAa— — = 2m(/x — h) — 4mAQ;-p = . 



Ai- 



A, 



(D9) 



a and A^/A^ are determined from Ea. (ID9[) . For simplicity, we consider the case of ^ — 0, for which one of the two 
SO split bands crosses the F point fc = 0. This is a typical situation which realizes the non-Abelian topological order 
(and hence the Majorana zero mode) as discussed in the previous sections. We obtain two sets of solutions of (|D9P : 
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(i) a = 7+, and A^/A-,- = (7_ - 7+)/(4mA). 

(ii) a = ^7-, and A^/A-^ = — 1(7- + 7+)/(4toA). 



Here 7± = \/ y/GAm^X^ + Am?h'^ ± Sm^A^. Thus, we have two solutions of the BdG equation for r < 



7+ 



4mA 



A^Zp{-f+r) 



Solution I 



(DIO) 



s'(^-l)''A^Z^_l(^7_r), = 



4mA 



(7+ +-/-)A^Zp{i-/-r). 



Solution II 



(Dll) 



Note that the solution I corresponds to the contribution from electrons with the finite Fermi momentum kp — 7+, 
while the solution II is dominated by electrons in the vicinity of the T point fc ~ 0. In the following, we try to construct 
a solution for the Majorana zero energy mode which is mainly formed by quasiparticles with k ^ corresponding to 
the solution II above. 

We now proceed to analyze the solution of the BdG equations (jDl|) and (ID2p for r > Vc- The dependence on 9 of 
the wave function can be separated by assuming the following form of the solution. 



^t(r,0) =e-*V>(r), u^(r, 0) = e^f />(r). 



(D12) 



Since Ua{T, ^) is multiplied by (—1)" when 9 is changed from to 2tt in this gauge^^l the solution (|D12p satisfies the 
correct boundary condition with respect to 6 only when the vorticity n is odd. Thus, we will find the zero energy 
solution only for odd n. 

Substituting ((DT2]) into Eqs.(lDl|) and (|D2|) . we obtain, 



J_ _^ ld_ 

2m \ dr"^ r dr 



{n - If 
4^2 



2A 



d n+1 



dr 



2r 



Af>* = 0, 



(D13) 



1/92 I g 

2m \ dr^ r dr 



{n + 1)^ 
4^2 



fi — h 



/i> + 2A(|:-\:^)/,> + A/r = o. 



(D14) 



In the following, we set ^ — 0, as mentioned before. We postulate that the solutions of (jD13|) and (|D14p consist of a 
slowly-varying function of r, g-\'{i){r) and the Bessel function; i.e. 



f^ir) = 5tM^^(aV), />(r) = gdr)Z,^i{a'r) 

with a' a constant. Substituting these expressions into (|D13p and (|D14p . we have, 

idZn^Qg Qg A *7* 7 ' 

2AZn+i — 2Aa Zri^gi — Aoi Z„+i + lig^Zn-i 

2 (jr 2 * * —5— 2 



,/2 



m dr dr 

1 dZ^ Qg^ 



2m 



- 2AZi^^ - 2Aa'Z2+i5t + ^glZl^ - hg^Z^^+i g^Z^^ = 0. 



(D15) 



(D16) 



(D17) 



m dr dr ^ dr j ^1 ■ -^i ^-t- 2 

To derive the third terms of the left-hand sides of Eqs. (jD16P and (jD17p . we have used the relations (jD6p . We assume 
that the Bessel functions Z^{a'r) appearing in the solution of /f^ are the first Hankel function H^\a'r). To solve 



Eqs. (jD16|) and (jD17P for g^^ and g^, we use the asymptotic form of the Hankel function _ffi (z) ^ ^;^exp[i( 
j{2y + 1))], and the asymptotic relations. 



1 dZJa'r) , V 

y I , . i ^ + -> 

Z^yar) dr r 



(D18) 



In the following analysis, it will be revealed that the parameter a' is determined as a pure imaginary number, (see 
Ea. (ID25p .) Thus, we approximate Z*^_i/Zn^ — > exp(i-|(n -I- 2)) in the asymptotic regime. Then, Eas. (ID16p and 

~2~ 2 

(ID17I) are rewritten into 



ia' dgf 
m dr 



2iX 



dgi 

dr 



2\a'g^ - Agle'^ 



hg^ 



2m 



5T = 0, 



(D19) 
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— ^ + 2^X^ - 2Xa'g, + Ag^e^f ("+^) - hg, - ^g, = 0. (D20) 
m ar dr ' 2m 

We here introduce new functions g±{r) = g^ ± ig^. From (jD19[) and (|D20[) . we have, 

- - 2^^"'5- + ^A5le^^ ("+^' + hg^ - ^g^ = 0, (D21) 

m ar dr Zm 

+ 2A^ + 2zAa'5+ - zA5le'^("+2) + /i<?+ - = 0. (D22) 

m dr dr 2m 

To solve Eas. (jD21[) and (|D22I) . we examine the fohowing two possible solutions; 

(a) g^ = 0, and g+ gives a nontrivial solution. 

(b) g^ = 0, and g^ gives a nontrivial solution. 

We, first, consider the solution (a). We postulate that 5+ — 5+ exp[i^(n— 1)] with 3+ a real function. Then, Eqs. (ID21[) 
and (|D22p are recast into 

d9+ . a' . , , 

- (D23) 

^5+ / ■ / ^-A 

These two equations are equivalent to each other provided that 

/i- A 



A 



(D25) 



Thus, we obtain a' as a pure imaginary, as already noticed above (below Eq. (jD18|) '). Therefore, the approximation 
used for the derivation of ((DT9| and from 0X6]) and ((DT7| is justified. Solving Eq.(|D23| or (|D24| with (|D25|) . 

we obtain the solution for 3+, 



.9+ 



= C+e/'^'^'-'T^. (D26) 



Then, from Egs. ljDlSp . and g^ — (?+/2, g^ — —igj^/2, we obtain the radial part of the wave functions and for 
the solution (a); 

/>(r) = ^e^K-De/'-'^^'^i/W, {^^-) > (D27) 



From the asymptotic form of the Hankel function H'^\iar) ^ y -^^^ exp(— ar — i-j{2v + 1)), we find that for large r, 
(|D27|) and (|D28p are given by 



P>, ^ .C+ / 2A fe-A. 



.C+ 2A 



h- A 



Thus, these functions are normalizable when /i — A > 0. We, now, consider the matching of the solution (a) for r > rc 
and the solution for r < rc at r = rc- As shown below, we can match the solution (a) for r > rc with the solution 
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II for r < Tc given by (jDll[) . Putting /3 = (n + l)/2 and Z^(i7_r) = H^\i^^r) in (jDll[) . we find the asymptotic 
behaviors of the radial part of the solution II for r < r^- 



f<{r) = A^h[% (z7_r) ~ A^e-^f("+i) 7^ 



(D31) 

7r7_r 



/<(r) = A^ijWi(*7_r) ^ -zA^e-'T(«+i) / e"^-"-. (D32) 

T- y 7r7_r 

To match the solutions, we adopt the following approximation. Assuming that the SO split is much larger than the 
Zceman energy, i.e. h <^ mX^, we have 7_ « h/{2\), A^/A^ ~ ^*(7+ + 7-)/(4wA) « — i. Then, noting that A — 
dX r ^ Tc, we can match the solution for r > Tc, (jD29p and (|D30p . with the solution for r < Vc, (|D31I) and (|D32p . by 
choosing C+ — 2-^2* exp(— + l))yl-f-. The solution for r < Vc, (|D3ip and (ID32I) . is not regular at r = 0, exhibiting 
logarithmic divergence ^ log(r) for r ^ 0. However, the solution is still normalizable. Thus, the solution (a) for 
r > Tc and the solution II for r < constitute the normalizable solution for the Majorana zero energy mode. This 
Majorana zero energy solution is constructed from quasiparticles in the vicinity of the V point fc '-^ 0, i.e. a single 
Dirac cone with a mass gap ~ /i, as mentioned before. 

We, now, examine the solution (b) for r > Tc] i.e. g+ = 0, and g- gives a nontrivial solution. In this case, we 
assume the solution of the form g-{r) = exp(i^(n + l))g_(r) with g-{r) a real function. Then, Eas. (jD2ip and (jD22[) 
are recast in 

^ = (D33) 
^^(— ' + ^).-. (D34) 

These two equations are equivalent when 



A 

For this choice of a', the radial part of the solution (b) is given by, 



a' ^i ^ ^ . (D35) 



/>(r) = ^e^ii^+^)e-rdr'^z^ ( t^-^r], (D36) 



We can not match (|D36[) and (|D37[) with the solution for r < rc at r ^ rc for any choice of the Bessel function 
Z,{i{A-h)r/X). 

Thus, when h > A and the vorticity n is odd, we obtain the only one normalizable zero energy solution which is 
given by (|DT2)) . ((D27)) . and ((D28)) for r > r^, and ((DTT|) for r < with /3 (n + l)/2 and = h!)^\ The field 
operator for the zero energy mode is given by 7^^ — J dr[u-i-{r)c!^{r) + uj^(r)cj^(r) + u|(r)c^(r) + £t|(r)cj^(r)], which is 
self-hermitian. Thus, the zero energy mode is a Majorana fermion. There is only one Majorana fermion mode in a 
vortex core with odd vorticity for h > A. 



Appendix E: Non-Abelian anyon in a time-reversal invariant s-wave superconducting state: Non-Abelian 

axion string'^*, and the Fu-Kane modeP^ 

In this paper, we mainly consider spin-singlet superconducting states under strong Zeeman magnetic field. Thus 
the time-reversal symmetry is explictly broken in the ground state. Indeed, the time-reversal breaking is necessary 
to obtain a non-zero TKNN number. From the bulk-edge correspondence, the non-zero TKNN number ensures the 
existence of topologically stable gapless Majorana fermions on a boundary and Majorana zero modes on a vortex. 
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However, it has been also known that non-Abehan anyons can be realized even when the ground state does not 
break the time-reversal invariance. Because of the time-reversal invariance, the TKNN number in this case is trivially 
zero. Nevertheless, the index theorem ensures the existence of a topologically stable Majorana zero mode in a vortex. 
This mechanism of non-Abelian anyon was discussed in refi^. Recently, Fu and Kane pointed out that Majorana 
fermions (and hence, non-Abelian anyons) realize in an interface between a topological insulator and an s-wave 
superconductors which preserves time-reversal symmetryj^^ The following analysis presents a general ground for the 
realization of non-Abelian anyons in such time-reversal invariant systems. 

To see this, consider the following Lagrangian for a 2+1 dimensional Majorana fermion TpM coupled with an scalar 
field (f> = $1 + $; 
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Here the Majorana fermion satisfies the Majorana condition, 

h2^M = V'M, 

and the Dirac gamma matrices 7^ and 75 are given by 



M- 



r 



cr^ 
at" 



75 = 




(El) 



(E2) 



(E3) 



where = (1, — tXi) and = (1, cr^) with the Pauli matrices cr^. We also suppose that $ is an s-wave condensate with 
the expectation value (<&) = $o- In sharp contrast to other non-Abelian topological phases, the system is time- reversal 
invariant, as was pointed out in rcf.^^. 

The above system has the following U{1) symmetry, 



M 



-20 



(E4) 



which is spontaneously broken by the condensate $0. Therefore, like an ordinary s-wave superconducting state, there 
exist a stable vortex solution that is given by 



^(x) = $o/(p)e'^ 



(E5) 



where p and (f) are the radial and angular coordinates from the vortex, respectively. The function f{p) vanishes on 
the core of the vortex and approaches to /(oo) = 1 far away from the core. 

The bound states in a vortex are studied by using the Hamiltonian of the system. 



n 



icjid,. 



(E6) 



and the unique zero mode, which satisfies "Huq — 0, is given by 



uo = C 



/ \ 

l + i 
1 - i 
\ / 



cxp 



-$o / drf{r) 




(E7) 



with a normalization constant The topological stability of the zero mode is ensured by the index theorem)^ 

From the Majorana condition (|E2I) in 2+1 dimensions, the operator of the zero mode 7 = / dxu\{x)il)yi{x) becomes 
real, i.e. 7^ = 7. Therefore, the vortex obeys the non-Abelian statistics. 

The above mechanism of non-Abelian anyons in an s-wave superconducting state is applicable to Axion strings 
in cosmological systems^i., and also to an interface between a topological insulator and an s-wave superconductor 
considered by Fu and Kane.— Indeed, identifying a gapless Dirac fermion on a surface of the topological insulator in 
the Nambu representation (-0^, ij)^^ "i/^j, — V'-f) and an s-wave Cooper pair due to the proximity effect with the Majorana 
field i/'M and the scalar field respectively, one can show that the BdG Hamiltonian considered in ref.— is essentially 
the same as the Hamiltonian ()E6|) . In this identification, the Dirac fermion in the Nambu representation satisfies the 
Majorana condition (|E2I) up to an unimportant factor. Furthermore, the electro-magnetic U{1) gauge symmetry in 
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the Fu-Kane model reduces to the U{1) axial symmetry (jE4[) . Therefore, for the same reason mentioned above, a 
vortex in the Fu-Kane model is found to obey the non-Abelian anyon statistics. 
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